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ABSTRACT 


This  study  is  concerned  with  the  mechanism  of  sand  movement  by  wind.  A 
method  for  calculating  the  rate  of  transport  was  developed.  This  method  was 
based  on  experience  gained  in  the  field  of  sediment  motion  in  rivers  and  all 
available  field  and  wind  tunnel  data  on  the  subject. 

It  has  been  found  that  the  basic  forces  causing  the  sediment  motion  are 
those  of  the  average  lift  L  and  the  fluctuating  part  L'  caused  by  the  tur¬ 
bulence. 

Another  factor  contributing  to  the  motion  in  the  case  of  wind  is  the 
effect  of  impact,  caused  by  the  particle  in  saltation,  in  disturbing  the  bed 
surface.  The  effect  of  impact  was  found  to  be  a  function  of  the  main  forces 
causing  the  motion  and  therefore  it  was  introduced  as  a  correction  for  the 
mean  lift  force  caused  by  the  distortion  of  the  fluid  field  around  the  bed 
particles. 

It  has  been  found  that  the  basic  principles  governing  the  rate  of  sediment 
transport  by  water  and  air  are  the  same.  The  only  difference  was  found  to  be 
the  effect  of  saltation  on  disturbing  the  bed  surface  in  the  case  of  air. 

The  results  of  this  study  are  represented  as  a  theoretical  relation  between 
the  flow  intensity  and  the  intensity  of  sediment  load.  The  effect  of  particle 
hidding  in  the  laminar  sub-layer  was  combined  with  the  impact  correction  to 
give  a  final  wind  correction  which  proved  to  be  a  function  of  the  parameter 
discribing  the  ratio  between  the  submerged  weight  of  the  particle  and  the  mean 
lift  force.  The  method  is  applicable  for  calculating  the  rate  of  sand  transport 
under  a  wide  range  of  wind  velocities  and  for  sand  sizes  ranging  from  a  0.145  mm 
to  1.00  mm.  Application  of  the  derived  method  for  calculating  sand  transport  by 
wind  from  natural  beaches  is  given. 
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The  problems  of  supply  and  loss  of  sediment  at  a  shore  line  are  of  con¬ 
siderable  importance  along  the  coast  line.  One  basic  mechanism  involved  in 

(25)* 

this  overall  problem  is  the  transportation  of  sand  by  wind  action  .  Many 
research  workers  have  studied  this  subject  in  the  laboratory  and/or  in  the 
field,  but  the  mechanism  of  the  sand  movement  by  wind  as  yet  has  not  been 
solved  completely.  One  of  the  principal  difficulties  is  that  the  forces  acting 
on  a  bed  particle  vary  with  respect  to  the  orientation  of  the  particle  in  the 
bed  and  the  different  stages  of  the  motion.  Particles  in  a  certain  size  in  a 
mixture  are  not  subjected  to  the  same  flow  velocities  as  they  are  in  the  case 
where  the  entire  bed  is  composed  of  material  of  its  own  size.  Also,  particles 
hidden  in  the  laminar  sublayer  are  not  affected  by  the  same  forces  as  those 
subjected  directly  to  the  main  turbulent  flow.  Most  of  these  difficulties 
were  overcome  in  the  case  of  sediment  motion  of  bottom  material  in  rivers 
after  extensive  studies.  In  1950  a  complete  theory  was  presented  by  Einstein^1^ 
which  permited  the  calculation  of  the  equilibrium  rate  at  which  various  dis¬ 
charges  will  transport  the  various  grain  sizes  of  the  bed  material  in  a  given 
channel.  Unfortunately,  in  moving  from  a  sand-water  system  to  a  sand-air 
system,  with  an  enormous  difference  in  density  between  air  and  water,  the  dis¬ 
turbance  of  the  bed  particles  caused  by  the  falling  grains  can  not  be  neg¬ 
lected  as  it  is  the  case  for  transport  by  flowing  water.  It  is,  however, 
believed  that  the  basic  principles  governing  the  motion  in  both  cases  should 
be  the  same. 

♦See  references  page  51. 
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With  the  recognition  of  the  effect  of  the  falling  grains  in  disturbing  the 
bed  particles  by  their  impact,  especially  at  low  transport  rates,  a  method  was 
developed  to  describe  the  mechanism  of  sand  movement  by  wind.  This  method  is 
based  on  the  Einstein  theory  for  sand  transport  by  flowing  water  and  all 
available  experimental  and  field  data  on  the  subject  of  sand  movement  by  wind. 

It  has  been  found  that  the  effect  of  surface  disturbance  caused  by  the 
falling  grains  can  be  introduced  as  a  correction  for  the  basic  forces  causing 
the  motion. 

It  has  been  also  found  that  if  the  Einstein  correction,  £ ,  for  particles 
hidden  between  larger  particles  or  in  the  laminar  sub-layer,  is  combined  with 
the  impact  correction,  a  final  correction  for  the  case  of  the  sand  transport 
by  wind  was  obtained.  This  combined  correction  proved  to  be  a  function  of 
the  parameter  \jr,  which  is  the  ratio  of  the  submerged  weight  of  the  particle 
and  the  mean  lift  force  caused  by  the  fluid. 

An  outline  of  using  the  proposed  method  for  calculating  the  transport  rate 
is  given.  The  application  of  the  method  was  checked  for  a  relatively  coarse 
sand  for  which  suspension  can  be  neglected. 

Summary  and  discussion  of  the  available  methods  for  calculating  the  rate 
of  sand  transport  by  wind  are  presented  next. 

II.  PREVIOUS  'WORK 

The  following  portion  of  this  study  deals  mainly  with  the  existing 
literature  on  the  subject  of  sand  transport  by  wind.  Information  on  the 
mechanism  of  transport  will  be  presented  as  given  by  the  various  authors. 
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Discussions  will  be  presented  in  the  following  part  of  this  study  ("statement 
of  the  problem") . 

(2) 

(i)  Types  of  sand  movement  by  wind  (after  Bagnold) 

For  values  of  U*,  the  shear  velocity  at  the  surface,  above  the  threshold 
value  U+t,  some  particles  from  the  bed  surface  are  put  into  motion.  The 

(2 

greater  the  value  of  U* ,  the  more  particles  will  be  put  in  motion.  Bagnold 
gave  three  possible  methods  of  movement:  saltation,  surface  creep  and  sus¬ 
pension. 

a)  Saltation 

(2) 

Bagnold  observations  ,  at  low  transport  rates  and  feeding  sand 
into  the  wind  tunnel  at  the  upwind  section,  show  that  the  main  motion 
of  the  grains  is  in  saltation.  Particles  rise  from  the  bed  with  neg¬ 
ligible  forward  speed,  are  accelerated  and  carried  forward  a  certain 
distance  by  the  fluid  flow,  and  finally  by  the  action  of  gravity,  they 
fall  back  to  the  bed  again. 

b)  Surface  creep 

Bagnold  described  the  mode  of  motion  by  surface  creep  as  follows: 

A  portion  of  the  energy  which  saltating  grains  have  gained  from  the 
wind  is  passed  on  to  the  grains  that  are  ejected  upward  to  continue 
the  saltation.  The  bulk  of  the  energy  is,  however,  dissipated  in  dis¬ 
turbing  a  large  number  of  surface  grains.  This  energy  is  ultimately 
all  lost  in  friction  between  the  surface  grains,  but  the  net  result 
of  the  continued  disturbance  of  the  surface  is  that  a  slow  forward 
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1.50  for  nearly  uniform  sand 

1.80  for  naturally  graded  sand 

2.80  for  sand  with  very  wide  range  of  grain  diameter 
The  general  characteristics  of  Bagnold  experiments  are  shown  in  Table  1. 

A  summary  of  his  measurements  using  0.25  mm  diameter  sand  are  shewn  in 


Table  2. 
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V  (3) 

b)  Kawamura 

Kawamura  following  the  same  basic  assumptions  as  Bagnold ,  but  using 
only  the  difference  between  the  shear  velocity  and  the  threshold  shear 
velocity  U*  ,  obtained  the  following  equation, 


q  =-  K, 


_f 

g 


(U*  +  U*t)4 


(U*-U*t) 


(2) 


where  K  is  a  constant  wH-ich  should  be  determined  experimentally  and  has 

XV 

been  found  by  Kawamura  to  be  2.78  for  a  sand  with  a  diameter  of  0.25  mm. 

Tables  1  and  2  show  the  summary  of  Kawamura  measurements. 

N  (5) 

c)  Zmgg 

Zingg  in  his  experimental  work  collected  sand  at  different  depths 
above  the  bed  for  a  known  period  of  time.  From  his  measurements  shown 
in  Tables  1  and  2,  for  five  different  grain  size  sand,  Zingg  obtained  the 


following  empirical  equation 

3/4 

q  =  C  <£■) 
z  Di 


r 


<V3/2 


(3) 


with  C  =0.83. 


z 

( 121 

d)  O’Brien  and  Rindlaubv 

These  investigators  performed  a  series  of  field  measurement  for  sand 
drift  by  wind  at  Clatsop  Beach  near  the  mouth  of  the  Columbia  River  in 
Oregon.  Wind  velocities  were  measured  at  elevations  between  0.25  and 
12  feet  above  the  ground.  The  rate  of  sand  movement  was  related  to  the 


wind  velocity  5  feet  above  the  sand  surface.  The  predominant  sand  size 


was  very  nearly  0.008  inch  (about  0.194  mm).  As  a  result  of  these  measure¬ 
ments  (Table  2),  they  obtained  the  following  formula. 

G  =  0.36  U53  (for  Ug  >  20  ft. /sec.)  (4) 

where  G  is  the  rate  of  movement  in  pounds  of  dry  sand  per  day  passing  an 
imaginary  line  1  ft.  in  length  drawn  perpendicular  to  the  wind. 

(4  6) 

e)  Experimental  studies  at  the  University  of  California,  Berkeley  ' 

Figure  1  shows  the  experimental  results  obtained  in  wind  tunnel  tests 
by  Bagnold,  Kawamura ,  Zingg  and  the  O’Brien  and  Rindlaub  field  measure¬ 
ments  upon  which  equations  1  through  4  were  formulated.  From  Figure  1 
it  is  clear  that  the  sand  transport  obtained  by  these  different  investi¬ 
gators  differs  widely  even  though  the  sand  considered  has  almost  the 
same  grain  size. 

In  order  to  reconcile  some  of  the  apparent  differences  in  the  various 

existing  relationships  for  the  rate  of  sand  movement,  studies  with  three 

different  grain  diameters  (Sands  A,  B  and  C)  were  conducted  at  the  Uni- 

(4) 

versity  of  California,  Berkeley.  The  characteristics  of  these  sands 
and  the  results  are  shown  in  Tables  1  and  2. 

From  the  results  of  these  studies  on  sand  transport  the  following 
conclusions  were  made:^^ 

1)  The  constant  in  equation  2  is  not  limited  in  range. 

2)  For  sand  B  it  was  impossible  to  find  a  value  for  these  constants 


which  would  permit  an  adequate  description  of  the  experimental  data. 


3)  The  Kawamura  equation  (equation  2)  includes  the  threshold  shear 
velocity,  U  ,  which  introduces  a  farther  uncertainty  in  the  cal- 
culations  of  transport  rates,  especially  since  this  is  influenced  by 

i)  personal  judgement  regarding  start  of  motion 

ii)  moisture  content  of  the  sand 

4)  The  O’Brien  and  Rindlaub  formula  (equation  4)  should  not  be  used 
for  calculations  of  transport  for  any  sand  diameter  except  of  the 
size  occurring  in  their  measurements. 

iii)  Shear  stresses  at  the  surface  (surface  drag) 

From  the  above  representations,  the  importance  of  the  shear  velocity  Ui 
at  the  bed  becomes  clear.  The  value  of  U  is  determined  by  measuring  the 
velocity  gradient  of  the  wind  near  the  bed.  For  a  turbulent,  steady,  uniform 
flow  of  fluid  over  a  stable,  rough  surface,  the  velocity  profile  can  be 
represented  by  the  prandtl  formula. 

0  *  c  Ioho  J-  (5> 

O 

where  C  is  the  rate  of  increase  of  velocity  with  log-height,  U  is  the  wind 
velocity  at  elevation  y  above  the  bed  and  y^  is  a  surface  parameter  character¬ 
izing  the  surface  roughness;  its  value  was  found  by  Bagnold  to  be  D/30  or 
0.033D. 

Under  steady  flow  conditions  over  flat  surfaces,  U+  which  is  equal  to 
the  square  root  of  the  shear  stress,^,  at  the  bed  divided  by  mass  density 
of  the  fluid,  pf,  has  a  very  important  physical  importance.  is  directly 


proportional  to  the  rate  of  wind  speed  with  log-height,  and  according  to  Von 

(16)  2.3 

Karman  ,  the  constant  of  proportionality  is  — —  ,  where  K  is  a  universal 

K 

constant  is  equal  to  0.4  for  conditions  of  no  sediment  motion.  Equation  5 
can  be  written  now  in  the  form 


U  =  5.75  D„  Log10  2- 

o 


(6) 


iv)  Velocity  distribution  over  a  drifting  sand  surface 

Once  the  wind  velocity  is  high  enough  to  move  sand  particles,  the  wind 

velocity  distribution  is  slightly  altered  by  the  sand  movement.  Plotted  on 

semi-log  paper,  the  velocity  distributions  remain  straight  lines,  but  as 
(2) 

shown  by  Bagnold  ,  they  all  seem  to  meet  at  a  certain  point,  which  he  calls 
a  "focus".  Bagnold  gave  the  velocity  distribution  above  a  surface  with  sand 
movement  by  the  equation 


U 


2.3 

K 


U  Log 

*  "^10 


\ 


u 


(7) 


where  uv  and  y'  are  the  coordinates  of  the  focal  point  for  drifting  surfaces. 

Zingg’s  experiments  with  five  different  sizes  of  sand^  showed  that  the 
projected  focal  points  (u1,  y')  appear  to  bear  a  relation  to  grain  size  in 
the  form, 


y'  =  10 

D 

(mm) 

in  mm. 

(8) 

and  u '  =  20 

D 

(mm) 

in  miles/hr. 

(9) 

The  range  of  grain  diameters  used  by  Zingg  ranged  from  0.20  mm  to  0,715 


9 


mm.  His  results  on  0.505  and  0.715  mm  sand  were  not  in  agreement  with 
equations  8  and  9. 

A  summary  of  all  available  measurements  on  the  coordinates  of  the  focal 

point  are  shown  in  Table  3. 

(5) 

Zingg's  measurements  showed  also  that  the  value  of  the  shear  stress 
at  the  bed  varies  from  the  value  obtained  indirectly  from  the  velocity  pro¬ 
file  with  K  =  0.4.  Zingg  proposed  a  value  of  0.375  for  K  in  equation  7, 
which  is  modified  now  to 

u  =  6.13  U+  Log10  -2-  +  u'  (10) 

III.  DISCUSSION  OF  THE  PREVIOUS  WORK  AND 
STATEMENT  OF  THE  PPQBLEM 

i)  General  remarks  on  Bagnold  and  Kawamura  Equations 

Although  both  the  Bagnold  and  Kawamura  equations  (equations  1  and  2)  are 

the  most  reliable  methods  at  the  present  time  for  calculating  the  rate  of 

sand  movement  by  wind,  they  are  based  on  assumptions  which  are  open  to 

question.  These  assumptions  can  be  summarized  as  follows: 

1  -  The  effect  of  turbulence  is  neglected  in  the  study  of  sand  movement  by 

wind.  Bagnold  mentioned  the  following, 

"The  effect  of  turbulence  is  not  appreciable  until 
far  higher  flow  values  are  reached" 

He  also  shows  indirectly  that  the  turbulence  is  greatly  responsible  for 
the  transportation  of  bed  material  by  the  measurements  shown  in  Figure  2-a 


- 
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and  describing  his  experimental  observations  in  a  special  wind  tunnel  as 
follows : 

"Sand  placed  at  the  mouth  itself  was  never  disturbed  even 
at  the  highest  wind  speeds  used,  despite  the  fact  that 
the  drag  and  normal  velocity  gradient  must  be  a  maximum 
here. " 

The  explanation  of  this  is  clear,  since  at  the  mouth  of  the  wind  tunnel,  tur¬ 
bulence  has  not  developed.  Figure  2-b  shows  the  sketch  for  the  boundary  layer 
and  the  turbulence  development  along  tie  wind  tunnel.  It  takes  a  distance  of 
about  4-5  mt.  (Fig.  2-b)  for  the  flow  to  be  fully  turbulent.  The  lack  of 
turbulence  at  the  mouth  of  the  wind  tunnel  (the  first  4-5  mt.)  can  be  the 
only  reason  for  the  failure  of  sand  to  move.  Since  Bagnold  has  indirectly 
proved  that  the  average  velocity  without  turbulence  is  definitely  not  able 
to  cause  the  movement,  it  must  be  concluded  that  the  turbulence  or  velocity 
pulsations  are  responsible  for  sand  movement  by  wind.  So  it  is  clear  now 
that  the  sediment  movement  by  wind  can  not  be  described  by  time  average 
values  alone,  and  more  understanding  of  the  motion  can  be  obtained  only  if 
turbulence  is  introduced. 

2  -  Both  Bagnold  and  Kawamura  considered  an  ideal  path  for  the  particle, 

namely,  the  saltation,  for  describing  the  rate  of  sand  transport.  A  variety 

of  arbitrary  assumptions  were  made  concerning  the  initial  and  final  velocities 

(2  3) 

of  the  saltating  grains.  ’  The  rise  of  the  grains  from  the  bed  was  assumed 
to  be  caused  by  the  mechanical  impact  of  the  falling  grains.  This  description 
applies  only  for  low  rates  of  transport.  However,  at  high  transport  rates 
the  whole  description  of  saltation,  for  initiating  the  motion  should  collapse 
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since  xt  is  only  the  direct  effect  of  air  which  puts  the  material  in  motion. 

(2) 

Bagnold  himself  mentioned  the  following  description: 

"At  low  wind  speeds  we  undoubtedly  have  saltation  and  the 
effect  of  grain  impact  predominating ,  at  higher  wind  speeds, 
with  an  appreciable  proportion  of  the  total  sand  flow  moving 
in  suspension  and  so  contributing  nothing  to  the  drag,  it 
is  possible  that  the  drag  due  to  the  remaining  saltation  is 
not  sufficient  to  keep  the  surface  wind  velocity  below  the 
threshold  value.  In  this  case  the  surface  wind  may  set 
grains  in  notion  by  its  direct  action  as  it  does  in  water." 

Therefore,  it  is  believed  that  the  description  of  the  transport  rate  based 
on  the  effect  of  saltation  only  is  misleading  and  should  be  explained  by 
the  variation  in  the  flow  conditions  at  the  bed.  I t  is  also  believed  that 
the  effect  of  saltation  especially  at  low  rates  of  transport  should  be  con¬ 
sidered  as  a  factor  contributing  to  the  initial  disturbance  of  the  bed  particle 
caused  mainly  by  the  turbulent  fluctuations  of  the  air. 

3  -  One  of  the  important  factors  which  should  be  introduced  in  the  study  of 
sediment  motion  is  the  grain  size  of  the  material.  Kawamura  ignored  this 
effect  completely  in  his  derivation,  assuming  that  the  same  equation  can  be 
used  for  any  grain  size.  Bagnold  assumed  that  tne  transport  rate  appears  to 
vary  approximately  with  the  square  root  of  the  grain  diameter.  Reduction  of 

(4) 

the  experimental  data  obtained  at  the  University  of  California  ,  including 
the  results  of  the  present  study,  are  plotted  in  Figure  3  which  shows  a 
family  of  curves  of  different  shear  velocities  for  the  change  of  the 
rate  of  transport  with  the  grain  diameter.  It  is  cleaa  from  Figure  3  that 
the  Bagnold  assumption  for  the  variation  of  the  transport  rate,  q,  with  the 
square  root  of  the  grain  diameter  (equation  1)  is  not  true.  Relations 
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obtained  in  Figure  3  are  of  great  importance  in  understanding  the  physical 
phenomena  of  sediment  movement  of  individual  grain  sizes  as  one  will  see 
later  in  this  study. 

4  -  All  the  available  equations  on  this  subject  contain  a  constant.  Sediment 
equations  which  depend  on  the  value  of  a  constant  are  not  desirable,  since 
this  constant  usually  does  not  have  a  physical  meaning.  Verification  of 
those  constants  proved  that  they  are  not  limited  in  range, 
ii)  The  Problem 

To  summarize  the  findings  of  the  above  discussions,  one  will  find  that: 

1  -  The  effect  of  turbulence  should  be  introduced  in  the  study  of  sand  move¬ 

ment  by  wind. 

2  -  The  effect  of  saltation  should  be  introduced  as  a  correction  for  the 

laws  governing  the  motion  and  not  as  the  main  cause  of  the  transport. 

3  -  The  effect  of  the  grain  size  of  the  bed  particles  should  be  considered 

as  an  important  factor  in  any  derivation  for  the  sediment  equations. 

4  -  The  basic  pictures  of  interaction  between  the  fluid  and  the  bed  particles 

was  not  discussed  by  any  investigator.  It  is  believed  that  the  problem 
can  only  be  solved  if  such  interaction  can  be  described  and  understood. 
Those  are  the  factors  upon  which  the  present  research  is  concentrated. 
The  problem  now  is  to  find  a  suitable  way  which  describes  the  overall 
mechanism  of  sand  movement  by  wind.  In  the  following  parts  of  this  report 
each  phase  of  the  phenomenon  will  be  explained  and  the  desired  relationships 


will  be  developed. 
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IV.  HYDRODYNAMIC  EFFECT  OF  THE  FLOW  ON  SEDIMENT  PARTICLES 


Numerous  studies  have  been  made  on  the  interaction  between  the  fluid 

flow  and  the  loose  material  forming  the  bed  of  a  stream.  The  results  of 

these  studies  led  to  the  development  of  theories  describing  the  phemenon  of 

sediment  transport  in  a  river.  One  of  the  most  reliable  methods  is  that  of 
(1  14  15) 

Einstein'  *  ’  .  Since,  in  both  cases  of  wind  and  water  transport,  one 

will  be  working  with  steady,  uniform,  turbulent  flow  passing  over  a  granular 
bed,  it  is  expected  that  a  basic  similarity  exists  between  the  motion  ol 
sediment  in  rivers  and  by  wind.  This,  of  course  does  not  imply  that  the  laws 
governing  the  same  phenomena  are  exactly  the  same,  nor  that  the  theories  de¬ 
scribing  the  motion  in  one  case  are  directly  applicable  to  the  other.  It 
is,  however,  reasonable  to  assume  that  some  fundamental  concepts  used  in  the 
derivation  of  one  theory  may  be  applied  in  the  derivation  of  the  other.  Since 
it  is  intended  to  use  in  the  present  study  some  of  the  basic  principles 
associated  with  the  motion  of  sediment  by  steady  streams  of  water,  one  could 
consider  it  appropriate  to  give  a  brief  outline  of  these  principles  and  how 
one  can  apply  them  to  the  present  problem, 
i)  Theories  of  sediment  transport  in  rivers 

In  a  study  of  the  history  of  developing  the  equations  and  methods  by 
which  sediment  transport  in  rivers  can  be  calculated,  one  will  find  that  it 

started  by  some  kind  of  equations  similar  to  those  now  in  use  for  sand 

(13) 

transport  by  wind.  For  example,  these  are  the  DuBoys  and  others  for  de¬ 
scribing  the  sediment  load  in  streams.  The  basic  concept  regarding  the 
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pattern  of  motion  was  that  the  loose  bed  is  sliding  in  layers  under  the 
action  of  the  flow  above.  No  effect  was  made  in  these  methods  to  explain 
the  actual  mechanism  of  interaction  between  the  solid  particles  and  the  flow 
field  General  information  therefore  can  not  be  deduced  from  these  methods 
and  their  applications  are  necessarily  limited  to  particular  conditions. 

Another  way  of  attacking  this  problem  was  to  consider  the  stability  of 

(23) 

the  individual  solid  particle  .  Unfortunately,  these  studies  did  not 
take  into  consideration  the  effect  of  turbulence.  As  a  result  the  flow  which 
is  responsible  for  the  forces  induced  on  the  particle  is  uniform  and  steady 
everywhere  on  the  bed,  causing  a  uniform  force  field.  Therefore,  a  particle 
that  starts  moving  at  some  point  of  the  bed  will  never  have  a  chance  to  come 
back  to  rest  at  some  other  point  on  the  bed,  a  mode  which  is  inconsistent 
with  the  actually  observed  mode  of  the  motion. 

(i  i4  is) 

In  the  last  twenty  years  a  new  method  was  developed  by  Einstein  ~  ’  ’ 

The  basic  concept  of  the  Einstein  theory  for  sediment  transport  in  rivers  is 
that  at  equilibrium  there  is  a  continuous  exchange,  at  the  same  rate,  between 
the  particle  in  motion  and  the  bed.  The  rate  of  deposition  is  found  to  be 
a  function  of  the  bed  load  rate,  while  the  rate  of  removal  of  grains  from 
the  bed  is  a  function  of  the  local  flow  intensity  and  the  probability  of  a 
particle  being  removed.  The  functional  relationship  between  the  "bed-load 
rate"  in  a  stream  and  the  flow  intensity  constitute  the  "bed-load  function” 
while  the  equation  expressing  this  relationship  is  called  the  "bed-load 
equation".  With  the  help  of  this  equation  it  is  possible  to  calculate  the 
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bed-load  rate  for  given  flow  conditions  and  bed  composition.  Einstein  extended 
his  work  on  bed-load  rate  a  step  farther  to  find  an  expression  for  the  sus¬ 
pended-load.  Since  true  suspension,  in  sand  transport  by  wind  is  impossible 
to  occur,  in  the  present  study  due  to  the  high  settling  velocities  of  sand 
particles,  information  on  suspended-load  is  not  needed.  However,  since  some 
of  the  basic  principles  used  by  Einstein  to  develop  his  bed  load-function  will 
be  used  in  the  present  study,  it  is  considered  desirable  to  give  a  comparison 
between  wind  and  water  transport, 
ii)  Comparison  between  air  and  water  transport 

From  the  above  general  description  of  the  Einstein  bed  load  function  one 
may  ask  the  question,  could  it  be  used  for  describing  the  motion  of  sand  under 
wind  action?  Before  attempting  to  answer  this  question,  one  should  consider 
the  principles  used  in  the  description  of  sediment  transport  which  led  to  the 
final  bed  load  equation  and  how  they  can  be  applied  to  the  present  problem. 

One  should  also  consider  the  enormous  difference  in  density  between  air  and 
water  in  making  the  comparison  between  the  eolian  and  aquatic  transport. 

Some  of  the  principles  and  assumptions  used  in  the  Einstein  theory  are*^ 

1)  The  turbulent  fluctuations  of  velocity  and  pressure  are  equally  important 
in  predicting  sediment  motion  as  the  average  value  of  the  main  flow. 

2)  The  "bed-load"  transport  directly  depends  on  the  grannular  material  and 
on  the  flow  pattei'n. 

(2> 

3)  The  motion  of  bed  particles  by  saltation  as  described  by  Bagnold  may¬ 
be  neglected  in  water. 
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4)  The  disturbance  of  the  bed  surface  by  moving  sediment  particles  may  be 
neglected  in  water. 

These  are  some  of  the  important  assumptions  used  by  Einstein  in  deriving  his 
bed-load  function.  The  most  important  of  them  is  the  first  one. 

The  important  role  played  by  turbulence  in  the  case  of  sediment  transport 
by  wind  has  been  demonstrated  in  different  parts  of  this  study.  Therefore, 
it  is  clear  now  that  the  principles  gained  in  this  regard  from  water-grain 
system  may  be  used  in  the  case  of  air-grain  system. 

Statements  2-4  inclusive,  require  comments,  since  movement  by  saltation 
and  its  effect  on  disturbing  the  bed  particles  and  the  transport  can  not  be 
neglected  in  the  case  of  air. 

Kalinske  has  shown  that  the  maximum  height  of  particle  "bounce"  is 

proportional  to  the  ratio  of  sand  density  to  the  fluid  density.  Therefore, 

for  equal  drag,  this  rise  of  grain  in  the  case  of  water  may  be  of  order  of 

— -  of  that  in  air.  Therefore,  neglecting  the  effect  of  saltation  in  the 
oUU 

case  of  water  is  justified,  since  the  particle  rise  will  be  in  the  order  of 
few  grain  diameters.  But  in  air,  on  the  other  hand,  particles  will  rise 
higher,  gain  momentum  from  the  moving  fluid  and  reach  the  bed  with  such  impact 
that  they  may  cause  considerable  disturbance  to  the  bed  surface. 

One  can  see  that  the  Einstein  description  of  motion  in  a  grain-water 
system  (namely,  "the  particle  moves  if  the  instantaneous  hydrodynamic  lift 
force  overcomes  the  particle  weight”)  should  be  re-examined  in  the  case  of 


grain-air  system  since  it  is  not  only  the  hydrodynamic  force  which  causes 
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the  particle  to  move  but  also  the  excitation  or  disturbance  of  the  bed  grains 
caused  by  the  falling  particles.  This  difference  should  be  expected  mainly 
at  low  transport  rates  where  the  effect  of  saltation  in  disturbing  the  bed 
particles  is  most  effective. 

The  effect  of  the  falling  grains  can  hardly  be  explained  by  a  simple 
equation,  since  it  depends  on  the  orientation  of  the  bed  particles,  the  rise 
of  the  falling  grains,  the  flow  intensity  and  many  other  factors  which  can 
not  be  enumerated.  Since  these  variables  can  not  be  measured  very  well,  the 
effect  of  falling  grains  must  be  determined  from  their  effect  on  the  motion  of 
the  bed  material.  The  significance  of  this  last  statement  will  become  clear 
from  the  following  derivations  of  the  sediment  transport  equation,  in  the 
case  of  wind,  based  on  the  Einstein  theory  for  sediment  transport  in  rivers, 
iii)  The  interaction  between  flow  and  bed  particles 

One  can  simply  look  at  the  general  problem  of  sediment  transport  as  a 
problem  of  a  steady  flow  over  a  granular  bed.  It  is  a  well  known  fact  that 
distortion  of  the  flow  field  around  a  solid  particle  resting  on  the  bed 
generates  a  lift  force  acting  on  the  particle,  even  if  the  latter  is  well- 
sheltered  within  the  sublayer.  This  force  will  tend  to  dislocate  the  particle 
and  move  it  away  from  the  solid  bed.  As  far  as  the  particle  is  still  m 
contact  with  the  bed,  the  lift  force  is  acting  only  vertically  upwards  and 
it  can  be  expressed  as 


L  =  C  o  —  a  n2 
L  wf  2  *  1 


(11) 


where,  C  is  the  lift  coefficient,  i.  is  the  instantaneous  velocity  acting  at 
L 

(17) 

a  distance  y  from  the  theoretical  bed.  Einstein  and  El-Samni  conducted 

flume  experiments  with  plastic  spherical  balls  0.225  feet  in  diameter  placed 

in  a  steady  stream  of  water.  The  theoretical  bed  has  been  determined  as  the 

reference  level  from  which  distance  should  be  taken  so  that  the  measured  values 

of  the  mean  velocity  would  give  the  best  fit  to  a  logarithmic  distribution. 

This  best  fit  was  obtained  when  the  theoretical  bed  was  taken  at  a  distance 

0.20  D  below  the  top  of  the  spherical  particle.  The  distance  Z  at  which  the 

o 

velocity  should  be  measured  in  calculating  the  lift  force  has  been  obtained 

simultaneously  with  the  determination  of  the  life  coefficient,  C  .  The  value 

Li 

of  the  lift  coefficient  was  obtained  by  measuring  the  pressure  at  two  places, 
one  at  the  top  of  the  spheres  and  the  other  near  their  bases.  The  pressure 
difference  which  is  a  measure  of  the  lift  force  can  be  expressed  as 

u2 

**  =  CL  Pf  T  (12) 

They  found  that  C  has  a  constant  value  of  0.178  for  a  wide  range  of  flow 

JL< 

conditions  provided  the  average  velocity  was  measured  at  a  distance  Z  = 

o 

0.35  D  from  the  theoretical  bed.  Since  in  both  cases  of  water  and  air  one 

will  be  dealing  with  a  steady  flow  of  fluid,  it  would  be  reasonable  to  assume 

that  the  value  of  C  is  about  the  same  for  air  and  water. 

L 

Unfortunately  the  problem  is  not  as  simple  as  mentioned  above.  Since, 
superimposed  over  the  above  conditions  one  has  the  effect  of  three  important 


factors , 


a)  Effect  of  turbulence 


b)  Effect  of  the  degree  of  sheltering  in  the  laminar 
sublayer 

c)  Effect  of  bed  particle  disturbance  caused  by  the 
falling  grains 

There  are  some  other  factors  but  of  less  degree  of  importance  than  the  above 
three.  Examples  of  these  are  the  particle  shape  ripple  formation  at  the 
surface  and  the  effect  of  the  sorting  coefficient  of  the  bed  material, 
a)  Effect  of  turbulence 

In  a  turbulent  flow  all  the  local  flow  parameters,  and  consequently  the 

local  lift  as  well,  vary  rapidly  with  time.  By  measuring  the  instantaneous 

values  of  the  lift  force  exerted  by  a  steady  stream  of  water  on  the  plastic 

spheres  mentioned  above,  El-Samni  found  that  the  total  lift  force  L  =  L  +  L’ 

behaves  like  a  random  variable  having  a  normal  distribution  with  mean  L  and 

s  tandard  deviation  a  =  L  n  =  — — —  . 

'o  2.0 

Knowing  the  distribution  of  L  +  L’,  one  can  establish  a  criterion  of 
stability  for  a  bed  particle  as  follows 

L  =  L  +  L  >  W,  (13) 

b 

where  is  the  submerged  weight  of  the  particle.  The  probability,  P,  that 


a  particle  resting  at  a 

certain 

location 

in 

the  bed  becomes 

just  ready 

move  can  be  obtained  as 

follows : 

P  =  P 

r 

[L  +  L 

>\  i 

(14) 

or  P  =  P 

r 

[•^T  > 
n  L 

‘o 

a* 
c r 

i 

1 

] 

(15) 
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Now,  using  equations  6  and  12  and  considering  the  velocity  0.35  D  above  the 

W 

bed,  as  found  bv  £l~Sammi^7^  the  ratio  — —  in  equation  15  can  be  written  as 

l  nQ 

»  <p„  -  p.)  g  D  A  ° 

_L.  =  - 5 - l - 1 - -  (16) 

L  %  %  '  CL  ’  Pf  [5'75  “»  1Og10  10'5}  L'  "  \ 


where  pg  and  pf  are  the  mass  densities  of  the  sand  and  air,  respectively. 


g  =  the  acceleration  due  to  gravity 


A^  =  shape  factor  for  the  grain  area 
=  shape  factor  for  the  grain  volume 


Equation  16  can  be  written  as 


Ps'pf  g  D 

- - -  - 7T 

pf  57 


no  CL  A1  5,75  (Log10  10 -5V 


Let  f  = 


Ps  -  Pf  g  D 


and  B+  = 


i)o  CL  Ai  (5. 75) 2  (Log10  10. 5)2 


and  since  -  has  a  normal  distribution  with  mean  zero  and  standard 

£'o 

deviation  0=1,  the  probability,  P,  in  equation  15  can  be  expressed  as 


p.  -i-  f" 

s**  b4  -  i 
% 

where  z  is  a  variable  of  integration. 
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b)  Hiding  effect 

Naturally  the  smaller  the  size  of  the  particles  D,  relative  to  the  thick¬ 
ness,  of  the  undisturbed  laminar  sublayer,  6,  the  less  pronounced  is  the  dis¬ 
tortion  of  the  fluid  field  around  the  particle,  and  the  intensity  of  uhe  lift 
force  will  be  reduced.  Also  when  one  considers  bed  particles  of  different 
grain  sizes,  it  can  be  realized  that  particles  of  a  certain  size  in  a  mixture 
are  not  subjected  to  the  same  flow  velocities  as  they  are  in  the  ca*e  where 
the  entire  bed  is  composed  of  material  of  its  own  size.  This  problem  was 
solved  by  Einstein^ 1  who  introduced  the  effect  of  particle  hiding  as  a 

correction  £  ,  for  the  lift  force  in  equation  13.  t  is  gxven  as  a  function 
D  , 

of  -  (Figure  4) ,  where  the  value  of  X  was  empirically  found  by  Einstein  to 

A 

be  1.39  5  for  a  smooth  bed  and  0.77  A  "for  rough  bed,  where  A  is  the  apparent 
roughness  parameter.  More  details  about  the  correction  £  is  presented  in 
reference  number  15, 

Now  equation  20  may  be  written  as: 

<v»  2 

i  r  -  ~ 

P  =  —  /  e  2.  dz  (21) 

42*  J  B*  *6  -  f- 
*  % 

c)  Disturbance  effect  caused  by  the  falling  grains 

Recalling  the  manner  by  which  the  Einstein  bed-load  equations  for  rivers 
was  derived,  it  can  be  seen  that  two  importance  assumptions  were  made,  namely, 

1.  The  motion  of  bed  particles  by  saltation  is  neglected. 

2.  The  disturbance  of  the  bed  surface  by  the  moving  sediment  particles 
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also  may  be  neglected  in  the  case  of  water. 

It  has  been  shown  m  part  (ii)  of  this  section  tnat  these  two  assumptions 
can  not  be  neglected  in  the  case  of  air  air-gram  system.  The  question  then 
comes  up  as  to  what  will  happen  to  the  lift  L  or  the  'i'  value  in  equation  21 
due  to  the  disturbance  action  caused  by  the  falling  grains?  Would  it  be 
necessary  to  introduce  a  correction  factor  for  that  effect?  The  answer  to 
these  questions  may  be  found  from  the  following  argument. 

If  one  considers  a  falling  grain  at  the  moment  it  makes  contact  with  the 
bed  particles,  it  can  be  visualised  that  the  grain  will  nave  an  appreciable 
velocity  at  the  moment  it  hits  the  bed.  Different  cases  can  be  considered 
as  follows: 

1.  The  grain  ricochets  upwards  from  a  solid  particle  supported  both  ways 

as  a  particle  in  position  (3)  in  figure  5.  It  will  rebound  with  a  velocity 

(2) 

nearly  equal  to  the  initial  velocity  of  impact  .  One  can  also  expect  tnat 
the  particle  in  position  (3)  will  become  more  exposed  tc  the  flow  than  before 
the  collision  with  particle  (5).  This  case  is  the  exceptional  one  except  for 
small  velocities  where  the  efficiency  of  the  collision  remains  hign. 

2.  The  falling  grain  will  help  put  other  grains,  not  completely  supported, 
in  motion  by  imparting  some  of  its  kinetic  energy  to  them  making  them  more 
exposed  to  the  flow.  Example'  of  these  are  particles  m  position  1,  2  and 

4  in  Figure  5.  In  the  majority  of  cases  the  impact  efficiency  will  be 

(19 ) 

relatively  low  and  the  gram  will  actually  rebound  to  only  a  small  height. 

Consider  now  particle  number  5  (Figure  5)  striding  the  particle  in  position 
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2,  for  example.  Particle  2  is  already  disturbed  by  the  direct  action  of  the 
wind.  It  is  very  difficult  to  say  how  the  impact  will  take  place  between 
particles  2  and  5,  but  generally  one  may  consider  the  particle  in  position  2 
and  try  to  determine  an  expression  for  the  effect  of  impact  caused  by  particle 
5  on  such  a  particle.  We  have, 

Force  =  rate  of  change  of  momentum 
(W  -  V  ) 

or  F  =  m  - ----- — —  (22) 

t 

where 

F  =  the  vertical  force  on  particle  2  as  the  result  of  a  collision  with 

falling  particle  5. 

m  =  mass  of  the  particle 

=  initial  velocity  of  particle  2:  assumed  =  o 

=  final  velocity  of  particle  2 

t  =  time  during  which  the  change  of  velocity  took  place. 

(2) 

W2  has  been  assumed  by  Bagiold  to  be  equal  to  B^  where  is  an  impact 

(3) 

coefficient  and  U  is  the  shear  velocity.  Kawamura  '  found  that  W  =  K 
*  *  2  1 

(U  +  U. , )  where  U.  .  is  the  threshold  shear  velocity  and  K,  is  a  constant. 
Thus,  it  is  reasonable  to  assume  that 

W2  =  C1  U,  (23) 

where  is  a  constant. 


t  can  be  looked  at  as  a  measure  of  the  time  required  for  the  replacement  of  a 
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particle  that  is  just  being  picked  up  by  the  flow  at  a  certain  spot  of  the 

bed  (particle  2,  Figure  5)  by  a  similar  particle  that  is  being  brought  to 

rest  at  the  same  spot  (particle  5,  Figure  5).  In  other  words  t  is  the 

Einstein  exchange  time^)  ,  i.e.,  t  can  be  assumed  to  be  proportional  to  the 

time  necessary  for  the  particle  to  settle  in  the  fluid  through  a  distance 

equal  to  its  own  size  in  the  fluid  at  rest.  If  the  settling  velocity  is 

denoted  by  V  then  one  can  write 
s 


t 


or  t 


K 


K 


JL 

V 

s 


gj 


JL 


2 


(24) 


(25) 


where  K  is  a  constant  of  proportionality  and  the  other  terms  are  as  defined 
before.  Now  substituting  equations  23  and  25  in  equation  22,  one  gets 


F  = 


C1  '  A2  D  V  X 


D  P. 


K 


2  '  pf  g 


) 


(26) 


P0  -  Pf  S  O  p  p  f 

or  F  =  K  [  -----  — r-  ]  [— - - -  l1'2  Fn  n  53 

O  P  o  tr  ^ 


-  ]  '  [U  D  2 


f  0. 


g  (Ps  -  Pf) 


(27) 


C1  A2 

where  K  =  — — —  =  constant 

3  K_ 


Using  equation  18,  equation  27  reduces  to 


P  1  /O  O  ^ 

F  -  K  *  [-2 - i - ]  [U  D  2  ] 

O  /  .  * 

g  (p  -  pf> 


(28) 
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Equation  (28)  shows  the  effect  of  surface  disturbance  F  to  be  a  function  of  the 
flow  intensity  ty,  as  given  by  equation  18  and  ( 


!>_-*_  ,i/*  oj  *  d3/2> 

P„  "  P,  * 


s  '  f 

Here  again  the  effect  of  the  force  F  on  the  sediment  motion  may  be  con¬ 
sidered  as  a  rapid  increase  in  the  effect  of  the  hydrodynamic  forces  acting  on 
the  bed  particles.  In  other  words  this  effect  may  be  considered  as  an  increase 

in  the  lift  coefficient  C  of  equation  11.  Therefore,  one  may  introduce  it 

L 

as  a  correction  for  the  lift  force  or  the  value  in  equation  21.  From  equation 

s  P  £  1/2 

28  one  would  expect  that  such  a  correction  should  be  a  function  of  ( - ) 

p  -  p. 

aJM.  s  f 


If  such  a  correction  is  called  I,  it  should  be  recognized  that  for  the 
case  of  a  water-grain  system,  the  value  of  I  will  be  equal  to  unity.  Since 
the  correction  I  is  of  opposite  effect  compared  with  that  of  the  correction  £ 
of  equation  21,  the  lift  force  should  be  multiplied  by  I  and  equation  21  may 
take  the  form 


P  = 


/ 


B*  *  ^ 


1 


(29) 


As  mentioned  before,  the  only  possible  way  of  describing  the  effect  of 
surface  disturbance  (correction  I)  is  by  its  effect  on  the  motion  of  bed 
material.  Therefore,  one  would  depend  on  experimental  measurements  to  find 
such  a  correction  as  it  is  seen  in  part  VII  of  this  study. 
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V.  THE  SEDIMENT-LOAD  EQUATIONS 


After  the  interaction  between  the  bed  particles  and  tne  fluid  flow  has 
been  explained  in  the  previous  part  it  remains  to  develop  a  method  which 
describes  the  final  solution  to  the  problem.  For  this  one  may  use  the 
equilibrium  condition  of  the  exchange  of  bed  particles  between  the  grains 
in  motion  and  the  bed.  For  each  unit  of  time  and  of  bed  area  the  same  amount 
of  sand  of  a  given  size  must  be  deposited  in  the  bed  as  are  scoured  from  it. 

Now,  consider  a  uniform  grain  size  sand  surface  of  unit  width,  over 
which  the  wind  is  blowing,  as  shown  in  Figure  6.  If  q  is  the  rate  of  sand 
moving  across  Section  (o)  in  dry  weight  of  sand  during  unit  time  and  unit 
length,  L  is  the  average  travel  distance  of  the  sand  particle  of  size  D  and 
t  is  the  time  consumed  for  replacing  a  particle  in  the  bed  by  a  similar 
particle  in  motion,  the  rate  in  dry  weight  of  sand  falling  on  a  unit  area 
of  sand  surface  per  unit  time  can  be  written  as 

°F  -  £  <30> 

The  number  N  of  particles  of  size  D  per  unit  area  of  the  bed  surface  that 
at  any  instance  becomes  free  to  move,  and  indeed  do  move  is  proportional  to 
the  probability  F  as  defined  by  equation  14  as  well  as  to  the  total  population 
of  similar  particles  per  unit  bed  surface,  therefore  N  can  be  written  as 


\  D 
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where  A^  is  a  shape  factor  for  the  particle  area.  The  rate  of  sand  scouring 
from  the  bed  in  dry  weight  per  unit  area  and  unit  time,  Gg ,  can  be  expressed  as 


O  A  D3 
s  g  2 


A2  P  Ps  8  D 
°r  Gs  *  t - 


where  A  is  a  shape  factor  for  particle  volume  and  the  other  terms  are  as 
defined  before. 

At  equilibrium  conditions,  G  should  be  equal  to  G  and  one  gets 

F  s 

A„  P  p  g  D 

l  =  ^  —~r —  (34) 

Einstein^  gives  the  average  travel  distance,  L,  by 


where  \  is  a  constant  and  the  other  terms  are  as  defined  above.  Now  sub¬ 
stituting  equations  35  and  24  in  equation  34,  equation  34  can  be  written  as 

P  A1  K2  ,  q  Pf  .1/2  1  1/2 

i-p  -  V  [  rr  ]  ,3W 

2  S  S  f  g  D 


and  A.  = 
* 


*=  — -  ( - - - )l/2 

V  g  ps  -  pf  g  D3 

A1  K2 
A2  X 


(38) 
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This  equation  is  identical  with  the  Einstein  equation  for  a  water-grain  system, 
except  the  correction  I  for  the  wind  condition  and  using  uniform  sand.  The 
practical  application  of  equation  (41)  necessitates  the  determination  of  the 
constants  A*  and  B+  and  the  |  and  I  values.  The  procedure  used  towards  this 
end  will  be  described  in  the  following  section. 

Now  if  one  defines  the  flow  intensity  as 


The  final  form  of  equation  (40)  will  be 


A*  ♦ 

1  +  A*  0 


dz 


(41) 


(42) 
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VI.  DETERMINATION  OF  A*  AND  B* 


The  application  of  equation  42  becomes  easy  once  the  values  of  A*  and 

B+  are  known.  One  method  of  obtaining  A*  and  B*  can  be  described  as  follows 

Suppose  that  a  rather  large  set  of  measurements  are  available  so  that 

another  set  of  values  as  defined  by  equation  18  and  <t>  values  as  defined 

by  equation  37  can  be  calculated.  Suppose  farther  that  those  measurements 

were  obtained  under  conditions  where  the  corrections  i  r 'd  I  of  equation  40 

are  known  to  be  unity  where  equation  40  takes  the  form 

2 


*  $ 


1  +  A*  ♦ 


^  / 

n~r  T' 


z 

2 


& 


dz 


(43) 


B*  * 


These  are  the  cases  under  which  the  basic  Einstein  theory  for  sediment  trans¬ 
port  was  derived. 

The  §  correction  is  known  to  be  unity  if  one  uses  in  the  above  set  of 
measurements  uniform  sediments  and  with  a  grain  size  coarse  enough  such  that 
D  >  >  1.  The  only  case  where  one  will  be  sure  that  the  correction  I  will  be 


h 


unity  is  the  case  of  a  water-grain  system  where  the  effect  of  the  disturbance 
caused  by  the  falling  grains  is  not  in  existence  as  described  before. 

Therefore,  one  car.  conclude  that  the  universal  cons  (.ants  A*  and  B+  can 
be  obtained  from  measurements  made  on  sand  transport  by  flowing  water.  As  a 
matter  of  fact  they  should  be  the  same  for  air  and  water  transport.  Einstein 
using  a  set  of  experimental  measurements  satisfying  the  conditions  mentioned 
above  obtained  the  following  values  for  A+  and  . 


BSSSSaaBE 


(i) 
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A*  = 


0 . 023 


(44) 


(29) 


B*  =  \ 

which  should  be  the  same  for  an  air-grain  system  as  explained  above.  The 


technique  used  for  finding  and  from  measurements  and  the  basic  theoretical 

( 1  24) 

equation  is  explained  elsewhere.  ’ 


VH.  DETERMINATION  OF  THE  CORRECTION  I 

As  mentioned  before,  the  best  method  of  finding  the  correction  I  is  from 
experimental  data.  Therefore  most  of  the  available  wind  tunnel  data  and 
field  measurements  on  sand  transport  by  wind  were  considered.  Table  1  shows 
the  main  characteristics  of  these  measurements  and  Table  2  shows  the  experi¬ 
mental  results  on  the  rate  of  sand  transport  q  and  the  shear  velocity  U*. 

All  these  data  were  used  in  finding  the  correction  I  except  those  of 
(3) 

Kawamura  .  The  Kawamura  experiments  were  conducted  in  a  wind  tunnel  of 
80  cm  height,  150  cm  length  and  5  cm  width.  Due  to  the  short  length  of  this 
tunnel,  it  is  clear  that  the  flow  did  not  have  a  chance  to  become  fully 
developed  and  the  shear  stress  2^  at  the  bottom  varied  considerably  along 
this  small  length.  This  fact  is  discussed  in  more  detail  in  text  books  on 
fluid  mechanics  where  the  development  of  the  boundary  layer  along  a  flat 
plate  is  discussed.  For  this  reason  the  Kawamura  data  were  excluded  from  this 
study. 

The  steps  necessary  for  obtaining  the  correction  I  are  as  follows- 
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1.  Equation  40  was  used  to  calculate  ©  for  different  —  values  with 

(17) 

q  =  0.5  as  found  by  El-Sarami  ,  =  43.5  and  =  0.1^3  as  found  by 

Einstein  and  using  the  tables  for  probability  integral.  Table  4  shows  these 
calculations.  Figure  7  shows  the  plotting  of  these  data. 

2.  The  value  as  defined  by  equation  18  was  calculated  for  all  the  available 
data,  the  shear  stress  U+  and  the  grain  diameter  were  obtained  from  Table  2. 

In  making  these  calculations  the  air  density  was  taken  as  1.22  x  10 

3  o 

gm/cm  (this  is  for  air  of  medium  humidity  at  15  C) ;  the  specific  gravity 

of  sand  was  assumed  to  be  2.65  for  all  sands  considered;  the  acceleration  of 

2 

gravity,  g  was  taken  as  980  cm/sec;  the  grain  diameter  D  was  considered  as 

D  Table  5  column  3  show  those  calculated  values  for  ty. 

5U 

3.  The  intensity  of  sediment  transport  0  as  defined  by  equation  37  was  cal¬ 
culated  using  the  available  experimental  values  on  the  rate  of  sand  transport 
q,  and  the  other  terms  in  equation  37  as  given  in  step  2,  The  calculated 
values  of  0  are  shown  in  Table  5,  column  4. 

4.  Now  using  the  theoretical  relationship  between  ©  and  ~  (Figure  7),  the 

\jr| 

values  of  —  can  be  obtained  for  each  corresponding  value  of  ©  obtained  from 

\lrf 

step  3.  Table  5  column  5  shows  the  values  of  —  as  obtained  from  Figure  7. 

5.  The  £  values  were  obtained  as  follows: 

a)  The  value  of  the  thickness  of  the  laminar  sub-layer  6  was  cal¬ 
culated  using  the  equation 


11.6  V 


5 


(45) 
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where  v  is  the  kinematic  viscosity  (taken  as  0.147  cm  /sec)  and  U+  is 
the  shear  velocity  which  can  be  obtained  from  Table  2. 

The  calculated  values  of  the  thickness  of  the  laminar  sub-layer,  5, 
are  shown  in  Table  5,  column  6. 

b)  The  values  of  were  calculated  as  shown  in  Table  5,  column  7. 

c 

c)  Since,  for  most  of  the  measurements,  the  bed  may  be  assumed  closer 

to  smooth  than  rough  (hydraulically) ,  the  values  of  X  =  1.39  5  can  be 

used  for  obtaining  the  I  value  as  recommended  by  Einstein^1^,  Therefore, 
°50 

the  values  of  — —  are  calculated  as  shown  in  Table  5,  column  8. 

A 

d)  Now  using  Figure  4  the  values  of  |  were  obtained  for  the  different 
D50 

values  obtained  from  above.  The  |  values  are  shown  in  column  9. 

Table  5. 

\!f  f 

6.  Finally  the  I  value  is  obtained  using  the  theoretical  value  of  —  =  , 

say,  obtained  from  step  4,  the  value  as  obtained  from  step  2  and  £  from 
step  5  using 


* 

I  =  ^|-  (by  definition) 


The  I  values  are  shown  in  Table  5,  column  10. 


As  was  mentioned  in  part  IV  of  this  study  and  from  equation  28,  the 

Pf  1/2  3  3/2- 

correction  I  appears  to  be  a  function  of  [ — — - — r  ]  [U+  D  ]. 

8  P-f* 

2  S  1 
Pg  P^  l/2 

Since  the  value  of  [ — - r  1  is  considered  a  constant  in  the  present 

*  <ps-p£> 

study,  the  values  of  I  as  obtained  above  and  the  corresponding  values  of 

3  3 /2  3  3 /2 

U*  D  were  plotted  in  Figure  8.  The  calculated  values  of  D 


for 


all  the  available  data  are  shown  in  column  11,  Table  5.  From  Figure  8  it  can 

3  3/2 

be  seen  that  an  inversely  linear  relation  exists  between  the  I  and  ij  D 
The  data  reasonably  fall  in  a  straight  line  with  a  slope  oi  minus  one  on 
a  log-log  plot  and  can  be  represented  by  the  equation 


3  3/2 

1  *  u*  D  =  a3 


where  A  is  a  constant. 

O 

This  relation  obtained  from  Figure  8,  can  be  regarded  as  one  of  the 
proofs  that  the  physical  picture  of  the  effect  of  the  surface  disturbance 
caused  by  the  falling  grains,  which  has  been  given  in  part  IV,  seems  to  be 
physically  sound. 


VIII.  COMBINING  OF  THE  EINSTEIN  CORRECTION  j  WITH 
THE  DISTURBANCE  CORRECTION  I 


In  using  the  Einstein  correction  g,  (Figure  4)  it  has  been  found  that 
most  of  the  available  measurements  on  sand  transport  by  wind  had  a  ^  value 

A 

less  than  unity.  Few  measurements  (sand  number  5  and  6,  Table  1)  had  a 
D 

value  of  -  >1.  For  those  few  measurements  it  was  found  that  if  the 

A 

Einstein  |  curve  is  modified  as  shown  in  Figure  4,  a  better  definition  of 


I  curve  (Figure  8)  was  obtained. 

As  previously  defined,  g  is  a  "hiding  factor"  which  expresses  the 
sheltering  effect  of  the  degree  of  submergence  in  the  laminar  sub-layer. 
Actually  a  change  in  the  ^  -  g  curve  can  mean  something  else  also.  The 
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nature  of  the  calculation  of  the  curve  j s  such  that  a hv  and  all  factors  not 
specifically  taken  into  account  in  some  other  manner  are  lumped  together  in 
the  one  i  value. 

Therefore  one  may  explain  the  modification  of  the  5  -  §  curve  in  the 

A 

present  study  by  some  factors,  not  yet  considered,  which  is  operative  in  the 
case  of  sand  transport  by  wind. 

Since  the  ~  -  §  curve  has  a  slope  of  2  on  a  log- log  plot,  it  can  be 

X. 

represented  mathematically  by  the  equation 


(|)2  ft  =  a4 


where  is  a  constant  and  X  =  1.39  5  as  given  before.  Now  substituting  the 

value  of  5  as  given  by  equa cion  45  and  considering  the  kinematic  viscosity 
v  as  a  constant;  equation  48  can  be  writ tc..  as 


3  2  2 
U*  D 

where  A,,  is  another  constant. 

5 

Using  equations  47  and  49  the  combined  correction  of  the  lift  force  in 

equation  (40)  can  oe  written  as 

3  3/2 

c  A  U*  D  ' 

I  =  -5.  ,  * -  (50) 

1  2  2 

U*  D  A„ 

*  3 


I  _  A  _JL. 

I  -  6  d1/2 


f 

1 


35 


~  Pf  1/2 

Since  (  -  )  was  considered  as  a  constant  in  the  present  study, 


equation  51  can  be  written  as 


5/1  = A/7  ■  \  t  ^  ru]1/2  ■  a7  *  ‘ 1/2 

S  JL 


where  A„  is  a  constant  and  as  defined  by  equation  18. 

t 

Equation  52  states  that  the  final  combined  correction  |  ,  which  may  be 
called  the  "wind  correction"  is  nothing  more  than  a  function  of  i{>'. 

Table  5,  column  12,  shows  the  j  values,  and  Figure  9  represents  the 
relationship  between  the  wind  correction,  j  ,  and  ty.  From  Figure  9  it  can 

t 

be  seen  that  a  linear  relation  exists  between  j  and  ijr.  The  data  reasonably 
fall  on  a  straight  line  with  a  slope  of  minus  0.6  on  log-log  paper  instead 
of  -0.5  as  indicated  by  equation  52.  The  difference  in  the  slope  can  be 
explained  by  the  scatter  of  the  data.  However,  one  can  fit  a  straight,  line 
with  slope-minus  0.5  through  the  data  of  Figure  9  with  a  reasonable  degree 
of  success. 

If  one  tries  now  to  summarize  how  the  transport  rate  can  be  calculated 
using  the  above  approach,  given  the  wind  speed  and  the  grain  size  of  the  sand, 
the  following  steps  may  be  used: 


1. 

2. 

3. 

4. 


Use  equation  7  to  find  the  shear  velocity  U, ,  i.e.,  U  =  6.13  U. 

y  y  ^  *. 

^10  r  +  u 


Calculate  the  value  of  using  equation  18,  i.e.,  t  = 
From  Figure  9  find  the  wind  correction  |  . 


Ps-Pf  g  D 


Knowing  the  value  of  t  j  find  from  Figure  7  and  then  the  rate  of 
sand  transport  q  using  equation  37,  i.e.,  0  =  q  .  Pf  1/2 ^  1  1/2 

v  %  \-pt  rs* 
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I X .  EXPERIMENTAL  VERIFICATION  USING  A 
RELATIVELY  COARSE  SAND 

1)  Scope  and  purpose 

The  experimental  program  consisted  of  two  series  of  tests  designated  as 
series  a  and  b.  The  main  purpose  of  these  tests  were  as  follows: 

a)  To  obtain  experimental  measurements  on  the  rate  of  sand  transport, 
q,  and  the  corresponding  value  of  the  shear  velocity,  IT  ,  using  a 
relatively  coarse  sediment  to  test  the  proposed  theory. 

b)  To  compare  the  experimental  values  of  the  threshold  shear  velocity 

0  ,  under  conditions  where  the  effect  of  surface  disturbance  is  elimi- 

*  U 

(9) 

nated,  with  the  corresponding  values  for  water  as  obtained  by  Shields. 

2)  Test  s-ries  (a) 

(i)  Experimental  Apparatus  and  Procedure 

The  experiments  were  conduced  in  a  wind  tunnel  located  in  Building  27 6 
at  the  Richmond  Field  Station  of  the  University  of  California.  The  tunnel  is 
4.0  ft.  wide,  1.60  ft.  high  and  100.0  ft.  long  and  is  constructed  of  plywood 
(Figures  10  and  11)  The  wind  was  generated  by  a  suction  fan  at  the  exit  end. 
The  mean  velocity  was  varied  from  20.0  to  70.0  ft. /sec.  by  controlling  the 
fan  speed. 

Wind  velocities  were  measured  using  a  standard  Prandtl  type  Pitot  tube 
which  was  attached  to  a  point  gage  and  introduced  into  the  air  stream  through 
the  top  of  the  wind  tunnel.  For  small  wind  speeds,  the  Pitot  tube  was 


connected  to  a  Magnehelic  gage  having  a  range  of  one-half  inch  of  water  and 


graduated  into  divisions  of  0.2  inch.  At  high  wind  speeds  the  Pitot  tube 
was  connected  to  an  Ellison  type  draft  gage  having  a  range  of  one  inch  of 
water  and  graduated  into  divisions  of  0.01  inch.  The  Magnehelic  gage  is 
usually  preferable  to  use  than  the  Ellison  type  gage  because  of  its  more 
rapid  response  to  pressure  change.  Figure  12-B  shows  both  the  Magnehelic 
and  Ellison  type  gages  used. 

Two  different  coarse  sands  were  used  in  this  series  of  experiments.  The 
mechanical  analyses  of  these  two  sands  are  shown  in  Figure  13  and  their 
characteristics  are  as  follows: 


Sand  Mean  grain  diameter  Sorting  Coeff. 

(mm) 


Grain  sizes  range 
(mm) 


D 


1.00 


1.20 


0.90  -  1.20 


E 


0.88 


1.41 


0.40  -  1.20 


The  sand  was  spread  over  a  length  of  62.00  feet  of  the  wind  tunnel,  with  a 
thickness  of  about  two  inches.  A  hopper  to  feed  sand  into  the  wind  tunnel 
automatically  was  placed  near  the  entrance  to  the  tunnel  (Figure  12-A) .  The 
rate  of  sand  feed  was  adjusted  to  be  equal  to  the  rate  of  sand  transport  as 
measured  by  the  sand  trap.  This  trap  was  8.0  ft.  long  and  consisted  of  18 
compartments  permanently  fixed  at  the  end  of  the  sand  bed.  In  order  to 
eliminate  the  side  wall  effects,  the  rate  of  sand  transport  was  measured 
only  over  a  width  of  2  ft.  in  the  central  part  of  the  wind  tunnel.  Sand  was 
removed  from  the  compartments  at  the  conclusion  of  each  run  with  a  vacuum 


cleaner.  The  time  consumed  in  each  run  and  the  weight  of  sand  collected  were 
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recorded.  After  each  run,  the  sand  surface  was  well  mixed  to  eliminate  any 
effect  of  sorting  and  was  leveled  to  be  ready  for  the  next  run. 
ii)  Experimental  Results  and  Discussions 

Velocity  distribution  on  the  sand  surface 

Vertical  wind  profiles  were  measured  at  a  distance  of  11  feet  upwind  from 
the  end  of  the  sand  bed  as  shown  in  Figure  11.  These  measurements  were  made  at 
the  center  of  the  wind  tunnel.  The  velocity  distributions  obtained  with  dif¬ 
ferent  fan  currents  are  shown  in  Figures  14  and  15  for  Sands  D  and  E,  respec¬ 
tively.  The  velocity  profiles  obey  the  logarithmic  formula  above  the  focal 
point  (Figures  14  and  15) .  The  focal  points  were  located  at: 

Sand  D  (D  =  1.00  mm) 

_ 50 _ 

y’  =  0.14  ft. 

U '  =  32.0  f t. /sec. 

These  do  not  agree  with  the  estimate  by  Zingg's  equation,  which  gives 
y*  =  0.0328  ft. 
u’  =  29.40  ft. /sec. 

Sand  E  (D  =  0.88  mm) 

_ 50  _ 

y'  =  0.125  ft. 
u'  =  32.0  ft. /sec. 

These  values  also  do  not  agree  with  Zingg's  estimate.  In  order  to  find 
a  better  definition  for  the  coordinates  of  the  focal  point  u'  and  y',  all 
the  available  experimental  data  shown  in  Table  3  and  the  data  obtained  above 
on  Sand  D  and  E  are  plotted  in  Figure  16.  Figure  16  shows  the  relation 


obtained  between  the  coordinate  of  the  focal  points  u'  in  ft. /sec,  y'  in  feet 
and  the  grain  diameter  in  mm.  Since  there  is  no  logical  method  for 
explaining  the  physical  meaning  of  the  focal  point  with  our  present  knowledge, 
it  is  proposed  that  the  .relations  obtained  in  Figure  16  may  be  used  to 
determine  u'  and  y'  in  equation  7  without  further  discussions. 

Rate  of  Sand  Transport 

The  amount  of  sand  caught  by  the  horizontal  trap  was  measured  for  shear- 
velocities  up  to  110  cm/sec.  The  shear  velocity,  UA ,  was  determined  by  the 
slope  of  the  velocity  distribution  in  Figures  14  and  15.  Table  6  shows  the 
results  of  this  series  of  experiments  using  Sands  D  and  E,  These  data  on 
the  shear  velocity,  UA ,  and  the  rate  of  transport,  q,  are  plotted  in  Figures 
17  and  18.  The  experimental  data  are  compared  with  Bagnold,  Kawamuxa ,  and 
the  proposed  method  as  follows: 

Sand  D  (D  =  1.00  mm) 

DU 

Figure  17  shows  a  comparison  between  the  experimental  data  and  the 
Bagnold  equation  (equation  1)  using  a  value  of  C  =  1.5.  No  value  for  the 
constant  in  the  Kawamura  equation  (equation  2)  was  found  to  represent  the 
experimental  data.  Table  7  shows  the  necessary  data  for  comparison  with  the 
proposed  method.  The  "If  values  were  calculated  from  equation  18,  the  ~  values 
were  obtained  from  Figure  9,  and  the  $  values  were  calculated  from  equation 
37.  Figure  9  shows  a  comparison  between  the  theoretical  and  experimental 
data,  the  results  are  in  good  agreement  if  one  considers  the  scatter  of 
data  shown  in  Figure  17.  So  it  can  be  concluded  that  the  proposed  method  can 
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be  used  for  a  relatively  uniform  and  coarse  sand  up  to  1.00  mm  diameter. 

Sand  E  (D,.,.  =  0.88  mm) 

_ 50 _ 

A®  shown  from  the  mechanical  analysis  in  Figure  13,  Sand  E  seems  to 
be  a  naturally  graded  sand  with  grain  sizes  ranging  from  0.40  mm  to  1.20  mm. 
Bagnold  recommended  a  value  of  1.8  for  the  constant  C  in  equation  1  for 
such  a  sand.  Figure  18  shows  a  comparison  between  the  experimental  data  and 
the  Bagnold  equation.  It  is  obvious  from  Figure  18  that  the  Bagnold  equation 
does  not  agree  with  the  experimental  values.  Here  again  no  values  for  the 
constant  K  in  Kawamura  equation  were  found  to  express  the  measured  data. 

K 

Einstein^  found  that  in  the  case  of  sediment  mixtures,  like  that  of 

Sand  E,  a  correction  factor  Y  should  be  introduced  to  describe  the  change  in 

K 

g 

the  lift  coefficient  in  a  mixture  and  is  a  function  of  —  as  shown  in  Figure 
21  (or,  of  the  Reynolds  number  of  the  flow  at  the  bed  surface).  The  length 
Kg  is  a  roughness  diameter  and  5  is  the  thickness  of  the  laminar  sub-layer. 

The  correction  factor  Y  is  unity  for  uniform  sediment.  Now,  introducing 
the  correction  Y,  the  flow  intensity  ty*  can  be  written  as 

'I'  =  —  Y  (53) 

*  i 

Table  8  shows  the  necessary  calculations  for  Sand  E  based  on  the  proposed  method 

and  introducing  the  correction  Y  as  recommended  by  Einstein^^  and  using 

K  =  D_„. 
s  oO 

Figure  9  shows  the  comparison  between  the  theoretical  and  the  experimental 


results.  It  can  be  seen  from  Figure  9  that  proposed  method  is  applicable  for 
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sand  with  a  wide  range  of  grain  sizes  if  the  correction  Y  is  introduced  and 
the  whole  calculation  is  based  on  the  mean  grain  diameter  of  the  sand  mixture, 


3)  Test  Series  (b) 

The  basic  difference  between  the  laws  governing  the  motion  in  the  case 

of  a  grain-water  system  and  a  grain-air  system  was  found  to  be  the  effect  of 

the  falling  grains  as  described  by  the  Correction  I.  If  the  above  statement 

is  true,  one  would  expect  that  the  laws  should  be  exactly  the  same  if  the 

effect  of  I  can  be  eliminated  in  the  case  of  sand  transport  by  wind.  The  only 

conditions  where  the  effect  of  the  falling  grains,  in  disturbing  the  bed 

particles,  can  be  eliminated  are  those  of  the  start  of  motion.  There  is  no 

(9) 

reason  to  doubt  that  the  general  principles  which  Shields  and  others  have 
applied  to  the  definition  of  the  condition  at  which  the  surface  grains  begin 
to  be  disturbed  by  the  fluid  flow  applies  for  all  grain-fluid  system.  But 
some  discrepancies  between  the  air  and  water  cases  are  rather  to  be  expected 
because : 

1.  Effect  of  falling  grains  (in  the  cases  of  feeding  s-  :d  at 
the  up-wind  section  of  the  tunnel)  is  not  the  same. 

2.  Exact  stage  at  which  initial  movement  may  be  said  to  have 
started  is  a  matter  of  personal  judgment. 

Shield's  experimental  values  of  the  dimensionless  function 

p  U  2 

a  =  [  — £ —  — 
p„  -  Pf  g  d 

z>  X 


(54) 
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in  terms  of  the  Reynolds  number  — —  for  grains  in  water  are  shown  in  Figure 

(2) 

20,  Bagnold  showed  that  the  corresponding  values  of  O.  found  for  air  gave 
a  curve  of  the  same  general  shape,  but  the  values  of  Ct  were  reduced  compared 
with  the  water  conditions.  In  Bagnold 's  experimental  determination  of  Or,  the 
effect  of  surface  disturbance  was  included.  It  is  believed  tnat  if  such 
effect  is  eliminated  the  values  of  or  for  water  and  air  cases  should  be  the 
same.  Therefore  the  purpose  of  this  investigation  is  to  find  the  values  of 
a  for  an  air-grain  system  under  conditions  where  the  surface  disturbance  can 
be  eliminated  and  to  compare  with  Shields'  curve  for  a  water-grain  system, 
i)  Experimental  apparatus  and  procedure 
Experiments  were  conducted  in  a  wind  tunnel  located  in  Building  160  at 
the  Richmond  Field  Station  of  the  University  of  California.  This  tunnel  is 
1  foot  wide,  60  feet  long,  and  1,28  feet  high  as  shown  in  Figure  22.  The 
channel  was  constructed  of  wood,  with  one  side  made  of  plate  glass  for 
observation  purposes.  The  wind  was  generated  by  a  blower  at  the  entrance  of 
the  wind  tunnel,  driven  by  an  A.C.  motor.  Wind  velocities  up  to  70  ft. /sec 
can  be  obtained  in  this  wind  tunnel.  Three  different  sands  were  used  Sand 
D,  E  and  Sand  B.  For  each  sand  the  test  procedure  was  as  follows: 

1.  The  first  15  feet  of  the  upwind  section  of  the  wind  tunnel  bed  was 
covered  by  a  fixed  artificial  roughness  (Sands  i),  E  and  B,  respectively). 

This  length  was  followed  by  one  square  foot  of  the  tunnel  of  loose  sand 
(Figure  22)  and  another  5  feet  of  fixed  roughness.  By  such  an  arrangement 
one  would  have  a  reasonably  fully  developed  turbulent  flow  over  the  loose 
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sand  and  due  to  the  small  area  of  the  loose  sand  the  effect  of  disturbance 
caused  by  any  moving  sand  can  be  eliminated. 

2.  The  wind  velocity  profile  at  which  motion  of  the  loose  sand  started 
was  recorded  and  the  value  of  the  threshold  shear  velocity  was  obtained. 

3.  The  value  of  CX  was  calculated  using  equation  54,  and  the  corresponding 
U 

*  D 

— —  was  calculated  for  each  sand. 

V 

iii)  Experimental  results  and  discussions 

The  results  of  this  series  of  experiments  can  be  summarized  as  follows: 


Sand 

(mm) 

U  (For  start  of  motion) 
(cm/sec) 

*D 

V 

a 

D 

1.00 

100.00 

68.00 

0.0462 

E 

0.88 

90.00 

54.00 

0.C43 

B 

0.30 

44.00 

9.00 

0.03 

These  data  are  plotted  in  Shields'  curve  (Figure  20)  -  where  the  agreement 
between  the  wind  experiments,  eliminating  the  disturbing  action  of  the  falling 
grains,  and  the  water  experiments  as  obtained  by  Shields  is  clear.  The  result 
of  this  series  of  experiments  proves  again  that  the  basic  laws  governing  the 
sand  motion  either  by  wind  or  by  water  is  the  same  and  the  only  difference  is 


the  effect  of  the  falling  grains  in  the  case  of  wind  condition. 
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X.  PRACTICAL  APPLICATION  OF  THE  METHOD 

This  part  further  explains  the  use  of  the  method  for  calculating  the 
rate  of  sand  movement  by  wind  and  the  use  of  the  necessary  graphs.  Sample 
calculation  for  sand  transport  from  a  beach  is  also  presented. 

a .  Application  to  uniform  sediment 

1.  The  wind  velocity  u  at  a  height  y  above  the  bed,  and  the  grain  diameter 
of  the  bed  particle  is  assumed  to  be  known. 

2.  Using  Figure  16  the  focal  point  u’  and  y'  can  be  obtained  and  from 
equation  10,  the  sfiear  velocity  U+  is  calculated. 

3.  The  parameter  is  calculated  from  equation  18. 

4.  Using  the  parameter  t  and  Figure  9,  the  wind  correction  j  is  obtained. 

5.  Now  with  the  know,  value  of  the  flow  intensity  the  <t>  value  is 

*  I 

obtained  from  Figure  7  and  hence  the  rate  of  transport  q  from  equation  37. 

b.  Application  to  sand  with  a  wide  range  of  grain  sizes 

The  same  procedure  as  for  uniform  sand  except  in  step  5  one  should  intro¬ 
duce  the  Einstein  correction  Y  which  can  be  obtained  from  Figure  21.  In  this 

case  the  flow  intensity  will  be  calculated  as 

* 

»*  -  *<f)  Y 

and  the  whole  calculations  will  be  based  on  the  average  grain  diameter  of  the 
bed  particles.  This  is  different  from  the  case  of  grain-water  system  where 
one  should  seperate  sand  into  seperate  size  ranges  and  calculate  them  in¬ 


dividually  . 
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c.  Sample  calculation  for  sand  transport  from  a  beach 

The  greatest  difficulty  in  applying  the  various  equations  and  graphs  for 
sand  transport  on  natural  beaches  is  the  basically  irregular  wind  conditions 
under  field  conditions.  Both  wind  duration  and  speed  changes  from  time  to 
time.  Each  reach  of  a  particular  beach  generally  is  different  from  every 
other  reach  (topography  and  grain  size  of  sand). 

Specific  procedures  and  calculations  for  determining  the  rate  of  sand 

(7) 

transport  inland  from  a  beach  by  wind  were  reported  by  the  writer  .  Those 
calculations,  using  the  Bagnold  equation  for  the  rate  of  the  sand  movement 
by  wind,  were  made  for  Salmon  Beach,  California. 

In  this  part,  application  of  the  outlined  method  will  be  made  for  reach 
No.  8  (see  reference  7)  of  Salmon  Beach.  California. 

Application  of  procedure  and  steps  of  calculation  as  applied  to  Reach  No.  8 
of  Salmon  Beach 

The  steps  in  calculating  sediment  movement  from  any  beach  using  the 
method  developed  in  the  above  discussion  may  be  illustrated  using  Reach  No.  8 
of  Salmon  Beach  as  follows: 

Step  1:  wind-speed  measurements  at  a  known  distance  above  the  bed  should  be 
obtained  and  the  direction  should  be  recorded.  Also  information  on  the 
dureation  of  this  particular  wind  speed  over  the  period  of  time  for  which  the 
amount  of  transport  is  to  be  calculated  should  be  recorded.  This  information 
was  obtained  for  Reach  No.  8  using  Reference  7.  For  Salmon  Beach  wind  speed 
measurements  were  made  18  feet  above  the  sand  surface  for  the  period  Sept.  1, 
1962  to  August  31,  1963.  Table  9  Column  (1)  show  the  wind  speed  as  recorded. 
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The  four  possible  directions  causing  inland  sand  transport  at  Salmon  Beach 

were  N,  NW,  W  and  SW  and  are  shown  in  Table  9.  Wind  duration  per  year,  t, 

and  length  of  reach  £  contributing  to  inland  transport  were  obtained  from 

r 

Reference  7  and  shown  in  Table  9. 

Step  2:  Grain  size  of  the  sand  bed  should  be  obtained.  For  Reach  No.  8  the 
mean  grain  diameter  is  0.58  mm. 

Step  3:  Using  and  Figure  16  the  coordinate  of  the  focal  point  for  the 

DU 

velocity  profiles  can  be  obtained.  For  Reach  No.  8  these  was  found  to  be 

y*  =  0.03  ft. 

u*  =18.00  ft/sec. 

Step  4:  For  each  wind  speed,  the  shear  velocity  U+  can  be  obtained  using 
equation  10;  namely 

U  =  6.13  U.  Log  —  +  u’ 

*  y 

The  U+  values  for  Reach  No.  8  are  shown  in  Table  9,  column  2,  using  the  above 
equation. 

Step  5:  Using  equation  18  the  value  can  be  calculated  using  U+  as  obtained 
from  step  4,  D  as  obtained  from  step  2  and  the  normal  values  of  p  ,  P  and  g. 

DU  S  X 

The  values  for  Reach  No.  8  are  shown  in  Table  9,  column  3. 

Step  6:  Using  Figure  9  and  the  value  as  obtained  from  step  5,  the  wind 

correction  |  can  be  obtained.  This  was  done  for  Reach  No.  8,  and  the  | 

values  are  shown  in  Table  9,  column  4. 

► 

Step  7:  The  flow  intensity  ^  =  '!>'  7  can  be  calculated  from  steps  5  and  6. 
.  *  X 
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Table  9,  column  (5),  shows  these  values  for  Reach  No.  8. 

It  is  clear  from  equation  42  that  for  values  of  ^  a  25  the  corresponding 
<J>  value  will  be  very  small  and  approaches  zero  i.e.,  there  is  no  motion; 
therefore,  one  would  suggest  that  the  calculation  should  start  from  higher 
wind  speeds  i.i  opposite  manner  of  that  used  in  Table  9  till  a  value  of 
=  25  is  reached.  Wind  speeds  corresponding  to  a  25  should  not  be  considered 
since  they  contribute  nothing  to  the  motion  of  the  bed  material.  It  is  believed 
that  t^is  is  a  very  good  characteristic  of  the  use  of  the  proposed  method  since 
it  tells  us  where  the  motion  starts  and  where  one  can  stop  his  calculation. 

This  fact  shows  how  the  Bagnold  formula  is  misleading,  since  it  gives  an  answer 
to  the  rate  of  transport  rate  for  any  value  of  U+  other  than  zero. 

Step  8:  Using  values  obtained  from  step  7  the  <P  values  can  be  obtained 
from  Figure  7.  Table  9,  Column  6  shows  those  values  for  Reach  No.  8. 


Step  9: 


Using  equation  37,  namely 


1/2 


1/2 

) 


The  rate  of  transport  q  can  be  calculated,  since  the  above  equation  can  be 


written  as 


q -♦*>.«>  t  ]1/2  (g  dV/2 

p£ 

or  for  Reach  No.  8  and  substituting  the  known  values  for  p  ,  p  ,  g  and  D  in 

f  s 


gm-cra-sec  system,  for  example, 


q  =  45.30  $ 


in  gm/cm  -  sec. 
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The  q  values  are  shown  in  Table  9,  column  7. 

Step  10:  Using  the  information  obtained  from  steps  1  and  9  the  total  transport 
rate  for  each  wind  direction  can  be  calculated  in  dry  weight  of  solids,  using 

Q  =  q  £  T 
r 

where  Z  is  the  length  contributing  to  transport  (step  1  and  reference  7)  in 
cm,  and  T  is  the  wind  duration  in  seconds  per  year.  The  values  of  Q  for 
Reach  No.  8  are  shown  in  Table  9,  column  9,  for  each  possible  wind  direction 
contributing  to  inland  sand  movement. 

Step  11:  The  total  annual  inland  transport  can  be  found  by  adding  the 

total  Q  for  each  wind  direction.  For  Reach  No.  8  of  Salmon  Beach,  this  was 

6  6 

found  to  be  2.36  16x10  lb/year  as  compared  with  5.40  x  10  using  the  Bagnold 
(7) 

formula  .  The  reason  for  this  big  difference  is  clear  from  the  discussion 
shown  in  step  7. 


XI.  SUMMARY  AND  CONCLUSIONS 

The  present  investigation  is  concerned  with  the  study  of  the  mechanism 
of  sand  movement  under  the  wind  action.  A  method  for  calculating  the  rate  of 
sand  transport  was  developed.  This  method  was  based  on  the  Einstein  theory 
for  sediment  transport  in  rivers  and  all  available  field  and  laboratory  measure¬ 
ments  on  the  subject  of  sand  movement  by  wind.  The  main  findings  of  this 
study  can  be  summarized  as  follows: 

1.  The  turbulence  or  velocity  pulsations  are  largely  responsible  for  sand 
movement  by  wind  and  their  effect  should  be  introduced. 


49 


2.  The  description  of  the  transport  rate  based  on  the  effect  of  saltation 
only  is  misleading  and  should  be  explained  by  the  variation  in  the  flow 
conditions  at  the  bed. 

3.  The  basic  forces  causing  the  sediment  motion  are  those  of  the  average  lift 
L  and  the  fluctuating  part  L'  caused  by  the  turbulence.  The  effect  of 
saltation,  in  the  case  of  wind,  can  be  introduced  as  a  correction  for  those 
basic  forces. 

4.  Basically  the  main  principles  governing  the  rate  of  sediment  transport  by 
water  and  air  are  the  same.  The  only  difference  appears  to  be  the  effect  of 
saltation  on  disturbing  the  bed  surface  in  the  case  of  air. 

5.  The  effect  of  surface  disturbance  was  introduced  as  a  correction,  I,  for 
the  main  lift  forces  causing  the  motion. 

6.  If  the  Einstein  correction  §  for  particle  "hiding"  in  the  laminar  sub¬ 
layer  is  combined  with  the  above  correction  for  surface  disturbance,  a  new 
"wind  correction"  was  obtained  which  seems  to  be  a  function  of  the  parameter 

given  by  equation  18. 

7.  The  conditions  of  start  of  sediment,  motion  as  given  by  Shields  in  the 
case  of  water  are  the  same  for  wind  conditions  if  the  effect  of  surface  dis¬ 
turbance  is  eliminated. 

8.  Experimental  curves  were  obtained  for  the  coordinate  of  the  focal  point 
u'  and  y  ’  as  a  function  of  the  grain  diameter. 

9.  The  method  is  applicable  for  calculating  the  rate  of  sand  transport  under 
a  wide  range  of  wind  velocities  and  for  sand  sizes  ranging  from  0.145  mm  to 


1.00  mm. 
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10.  The  derivation  of  the  method  is  based  on  uniform  sand.  It  has  been  found 
that  it  can  be  used  for  sand  with  a  wide  range  of  grain  sizes  if  the  Einstein 
correction  Y  is  introduced.  Practical  application  of  the  method  for  uniform 
and  non-uniform  sand  is  given  together  with  a  sample  calculation  for  sand 
transport  from  a  beach  by  wind  using  the  derived  method. 
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XI I II.  FIGURES 


RATIO  OF  SAND  TRANSPORT  (  q  in  gram /cm- 


I 


(o)mean  air  speed  430  cm/sec  (below  threshold)  V*=i95  ;  (b)  mean  air  speed 
490  cm /sec,  V*=36  ;  (c)  mean  oir  speed  930  cm/sec,  V*.=  92. 
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Fig  2 a  variation  of  the  flow  of  sond  with  the  distance  along 

the  tunnel  after  Bagnold 
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Fig.  2b  development  of  turbulence  and  the  boundry  layer  along 

the  wind  tunnel 
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FIG.  5  EFFECT  OF  FALLING  PARTICLE  (5)  IN  DISTURBING 
THE  BED  PARTICLES 


FIG.  6  EXCHANGE  OF  BED  PARTICLES  WITH  PARTICLES 

IN  MOTION 
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FIG.  12-A  AIR  INTAKE  ft  HOPPER  FOR  SAND  FEEDING 


FIG. I  2-  B  MAGNETIC  ELLISON  GAGES  AND  PITOT 
TUBE  FOR  VELOCITY  MEASUREMENTS 
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FIG.  13  MECHANICAL  ANALYSIS  OF  SAND  D  8  E 
(  PRESENT  STUDY) 
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FIG.  14  VELOCITY  DISTRIBUTION  ABOVE  SAND 

SURFACE  WITH  SAND  MOVEMENT  {  SAND  D  ) 
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FIG.  15  VELOCITY  DISTRIBUTION  ABOVE  SAND  BED 
(WITH  SAND  MOVEMENT)  SAND  E 


FIG.  16  RELATIONSHIP  BETWEEN  U  AND  y  TO  AVERAGE  GRAIN  DIAMETER 
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FiG.  17  COMPARISON  BETWEEN  EXPERIMENTAL  RESULTS  AND 
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FIG.  18  COMPARISON  BETWEEN  EXPERIMENTAL  RESULTS  AND 
BAGNOLD  FORMULA  (SAND  E) 


FIG.  19  COMPARISON  BETWEEN  EXPERIMENTAL  RESULTS  AND  PROPOSED  METHOD  (SAND  D,E) 
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FIG.  21  CORRECTION  Y  FOR  CHANGE  IN  THE  LIFT  COEFFICIENT 
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FIG.  22  WIND  TUNNEL  USED  FOR  SERIES  b  AND  LAYOUT  OF 
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Table  1.  Characteristics  of  the  Available 
Experimental  and  Field  Data  on  Sand 
Transport  by  Wind. 

„  wind  tunnel  dims. 


Sand 

Tested  by 

date 

Sorting 
Ccsff .SQ 

°50 

(mm) 

Length 

(feet) 

Width 

(feet) 

Height 

(feet) 

1 

Bagnold 

1936 

1.00 

0.250 

30.0 

1.00 

2.50 

ft2 

2 

Zingg 

1942 

- 

0.200 

56.00 

area  =  3 

3 

Zingg 

1942 

- 

0.2.75 

56.00 

area  =  3 

ft2 

4 

Zingg 

1942 

- 

0.360 

56.00 

area  -  3 

ft2 

5 

Zingg 

1942 

- 

0.505 

56.00 

area  =  3 

ft2 

6 

Zingg 

1942 

- 

0.715 

56.00 

area  =  3 

ft2 

7 

Kawamura 

1951 

1.00 

0.250 

4.50 

0.165 

2.60 

8 

O'Brien  & 

1935 

- 

0.20 

field 

measurements. 

9 

Rindlaub 

Horikawa 

1960 

1.15 

0.20 

60 

1.00 

1.28 

A 

Bel ley 

1962 

1.23 

0.44 

100 

4.0 

2.50 

B 

Belley 

1962 

1.15 

0.30 

100 

40 

2.50 

C 

Kadib 

1963 

1.24 

0.145 

100 

40 

2.5 
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Experimental  and  Field  Data  on  the 
Sand  Transport  rate  q  and  the 
Shear  Velocity  U^. 


Sand  number  1 

Bagnold,  1936 

D  =  0.25  mm 
50 


Sand  number  2 

Zingg  1941 

D_-  =  0.20  mm 
50 


Sand  number  3 

Zingg  1941 

D_rt  =  0.275  mm 
50 


u* 

q 

u* 

q 

u* 

q 

cm/sec 

gm/cm-sec 

cm/sec 

gm/cm-sec 

cm/sec 

gm/cm-sec 

42.0 

0.0382 

60.00 

0.187 

19.2 

0 

46.3 

0.0895 

70.00 

0.304 

25.0 

0.029 

70.50 

0.217 

86.50 

0.505 

40.40 

0.118 

86.5 

0.603 

89.50 

0.83 

50.50 

0.25 

52.6 

0.131 

42.50 

0.075 

62.00 

0.44 

65.00 

0.262 

58.00 

0.167 

88.00 

1.22 

71.50 

0.217 

63.60 

0.30 

73.00 

0.342 

74.50 

0.605 

79.50 

0.354 

85.00 

0.78 

80.60 

0.572 

45.60 

0.054 

Sand  N.4 
Zingg  1941, 

D_  =0.36  mm 

50 

Sand  N.5 
Zingg  1941, 

D„  =0.505  mm 

50 

Sand  N. 
Zingg, 

6 

1941  D  =0.715  mm 
5U 

u* 

q 

U* 

q 

U* 

q 

cm/sec 

qm/cm-sec 

cm/sec 

qm/cm-sec 

cm/sec 

qm/cm-sec 

58.0 

0.284 

50.2 

0.187 

67.0 

0.545 

81.00 

1.27 

62.00 

0.329 

80.0 

0.746 

102.0 

2.28 

63.50 

0.42 

83.0 

1.03 

77.50 

0.975 

73.00 

0.67 

101.0 

1.67 

97.00 

1.92 

81.00 

0.905 

62.0 

0.48 

77.00 

0.738 

76.5 

0.610 

74.50 

0.575 

60.5 

0.319 

68.00 

0.475 

91.0 

1.29 

60.00 

0.385 

* 


[; 
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Trt}  i q  o  Oon tinusd 


Sand  N. 7 

S 

* 

and  8 

Sand  9 

Kawamura 

1951 

0  'Brien 

and  Rindlaub, 

Horikawa , 

i960 

D_  =  0. 
50 

25  mm 

D50= 

0.20  mm 

°50  ‘  °'20 

mm 

U* 

q 

U* 

q 

q 

cm/sec 

gm/cm-sec 

cm/sec 

gm/cm-sec 

cm/sec 

gm/cm-sec 

27.0 

0.01 

37.00 

0.378 

26.8 

0.0006 

33.0 

0.13 

37.20 

0.278 

40.00 

0.064 

36.0 

O.oO 

36.20 

0.344 

50.00 

0. 140 

38.0 

0.22 

34.50 

0.212 

60.00 

0.23 

42.0 

0.24 

25.00 

0.0875 

70.00 

0.33 

46.0 

0.40 

23.80 

0.074 

80.00 

0.52 

54.0 

0.65 

21.50 

0.0541 

90.00 

0.64 

65.0 

1.15 

20.6 

0.0464 

100.00 

0.83 

73.0 

1.60 

81.5 

2.20 

87.0 

2.75 

97.0 

3.45 

109.0 

4.30 

Sand  A 

D  =  0 

50 

Bel ley  1962 
,44  mm 

Sand  B, 
°50  =  ° 

Belley  1962 
.30  mm 

Sand  C, 

D  =  0. 

50 

Kadib  1963 
145  mm 

u* 

q 

q 

U* 

q 

cm/sec 

gm/cm-sec 

cm/sec 

gm/cm-sec 

cm/sec 

gm/cm-sec 

22.0 

0 

24.40 

0.0017 

27.80 

0.0277 

33.60 

0.059 

30.0 

0.012 

16.00 

0 

36.80 

0.088 

35.0 

0.105 

33.00 

0.25 

41.0 

0.144 

38.0 

0.182 

38.00 

0.32 

59.5 

0.189 

39.0 

0.220 

41.00 

0.36 

62.6 

0.270 

41.0 

0.232 

44.00 

0.39 

64.50 

0.365 

46.0 

0.380 

49.00 

0.50 

73.00 

0.475 

50.0 

0.506 

59.0 

0.74 

76.00 

0.53 

55.00 

0.780 

66.0 

0.93 

79.00 

0.649 

64.00 

1.180 

80.50 

0.736 

84.00 

0.95 

♦These 

data  are  field 

measurements 

,  U*  was  calculated 

using  equation  7  , 

since  q 

was  given  as 

a  function  of 

wind  speed  5  ft.  above  the 

bed. 
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To  K1  r»  Q  Cl » mm  o  wt  o-P  4-V»  /\  At»nn  1  aUI  a  *ma«%  1 

AM  VAV  V»  »  W  MI1UHM  A  j  va  witv  n  VOdkXHUXC  UAJJCJ.  X1UC11  UCXJL. 

data  on  the  Coordinates  of  the  Focal  Point  u'  and  z'. 


Tested 

D50 

U' 

3' 

Sand 

by 

(nun) 

(ft/sec) 

(ft.) 

2 

Zingg 

0.20 

5.90 

0.0075 

3 

Zingg 

0.275 

7.35 

0.010 

4 

Zingg 

0.36 

11.00 

0.012 

5 

Zingg 

0.505 

14.60 

0.02 

6 

Zingg 

0.715 

26.00 

0.070 

A 

Bel ley 

0.440 

13.00 

0.0144 

B 

Belley 

0.30 

9.00 

0.010 

C 

Xadib 

0.145 

6.40 

0.0125 

SE^. 


81 


\ 


lire 

Table  4.  Calculations  for  —  =  vs.  <5> 
relation  as  defined  by  equation  (40). 


ilrf 

Using  equation  40  one  can  calculate  <J>  for  chosen  —  values  and  using  the  table 
of  probability  integrals  using 


A.  =43.5 
¥ 


B.  =  0.143 
* 


and 


2.00 


*  -« 

*  I 

*  B  -  — 

*  *  n 
‘o 

n  -z2 

/  oo  e  dz 

J  i 
♦,b«-  r- 

‘o 

20.0 

2.86 

0.86 

0.1952 

0.0055 

15.00 

2.14 

0.14 

0.4442 

0.0186 

10.00 

1.43 

-0.57 

0.716 

0.058 

5.00 

0.715 

-1.285 

0.903 

0.24 

1.00 

0.143 

-1.857 

0.968 

6.95 

Table  5.  Necessary  Calculations  for  the  Corrections  £ ,  I  and 
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(present  study) 


Sand  D  (D^ 

=  1.00  mm) 

Sand  E  (D,.. 

50 

=  0.88  mm) 

U* 

q 

U* 

q 

cm  sec 

gin  cm-sec 

cm/sec 

gm/cm-see 

28.5 

0.46 

48.75 

0.215 

37.00 

1.71 

71.60 

1.31 

37.50 

1.00 

60.00 

1.36 

40.0 

1.40 

78.60 

2.37 

45.0 

1.65 

100.00 

4.20 

51.00 

2.60 

107.0 

4.75 

61.00 

4.10 

63.6 

0.59 

62.00 

4.00 

70.0 

1.925 

64.00 

4.50 

76.20 

2.87 

72.00 

5. 15 

85.0 


6.74 
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npfprmi  Tin  nil  fif  A  f  vnm 


1 3 O Xfcf  I  •  yc  iciuuua  tluu  v/x 


T* 


Experimental  Data 
(Sand  D) 


Sand  D 

D_rt  =  1.00  mm 
50 


ps  "  pf 


gD  =  2.14X105,  1/  =  2-‘iy1^- 
3/2 


d3/2  =  =  0.0316  cm3/2 

*  -  -375  259  *  1<r4 


=  8.20  x  10  3  q 


u* 

cm/sec 

q 

gm/cm-sec 

* 

i* 

I 

*  f  =  % 

<t> 

48.75 

0.215 

90.10 

0.25 

22.50 

1 . 76xl0~3 

71.60 

1.310 

41.80 

0.42 

17.50 

1. 07x10 ~2 

60.00 

1.360 

59.40 

0.34 

20.5 

1 . 12xl0~2 

78.60 

2.370 

34.60 

0.49 

16.90 

1 . 95xl0~2 

100.00 

4.200 

21.40 

0.68 

14.50 

3 . 45xl0_2 

107.00 

4.750 

18.20 

0.75 

13.60 

3 . 9xl0~2 

63.60 

0.590 

52.9 

0.365 

19.3 

4.84xl0”3 

70.00 

1.925 

43.6 

0.410 

17.9 

1 . 60xl0~2 

76.20 

2.87 

36.9 

0.47 

17.3 

2. 36x10 “2 

*  The  values 

-  t 

of  -  were 

obtained  from  Figure  9. 
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Table  8.  Determination  of  ^  and  <P  to  be 
Compared  with  the  Theoretical  Values 
(Sand  E) 

D__  =  0.88  mm 
50 


u* ** 

cm/sec 

q 

qm/cm  sec 

* 

.  * 
i/I 

6 

(mm) 

D50 

& 

** 

Y 

X 

0 

28.5 

0.46 

230.0 

0.132 

0.600 

1.47 

0.74 

22.00 

4 . 52xl0~3 

37.00 

1.71 

137.0 

0.180 

0.462 

1.91 

0.62 

15.20 

1.685x10  ' 

37.50 

1.00 

134.0 

0.190 

0.455 

1.94 

0.62 

15.80 

9.86x10  3 

40.00 

1.40 

117.50 

0.210 

0.428 

2.06 

0.60 

14.90 

1.38x10  2 

45.00 

1.65 

92.50 

0.255 

0.380 

2.32 

0.57 

13.50 

1.625x10 

51.00 

2.60 

72.50 

0.300 

0.336 

2.63 

0.55 

12.10 

2.56x10  2 

61.00 

4.10 

50.6 

0.380 

0.280 

3.15 

0.54 

10.60 

4.08x10  2 

62.00 

4.00 

49.00 

0.390 

0.276 

3.18 

0.54 

10.30 

3,98x10  2 

64.00 

4.50 

46.00 

0.410 

0.267 

3.30 

0.54 

10.20 

4.49x10  “ 

72.00 

5.15 

36.3 

0.48 

0.238 

3.70 

0.53 

9.25 

5.14x10  2 

85.00 

6.75 

26.0 

0.60 

0.201 

4.38 

0.53 

8.30 

6.65x10  2 

*  The  values  of  -  were  obtained  from  Figure  9 

**  The  Y  values  were  obtained  from  Figure  2i 
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Table  9.  Calculations  for  Sand  Transport  (in  land) 
using  the  Proposed  Method  for  Reach  8,  Salmon  Beach,  California 


U18  ft. 

U* 

( h 

q 

ft/sec 

cm/sec 

to 

(3 

4> 

gm/cm-sec 

(1) 

(2) 

(3) 

(6) 

(7) 

18.2 

0.40* 

>25 

0 

0 

18.8 

1.60* 

>25 

0 

0 

20.00 

4-00* 

0> 

u 

>25 

0 

0 

21.00 

6.00* 

a  a) 
o 

>25 

0 

0 

21.80 

7-60* 

G  G 
O  H 

>25 

0 

0 

22.90 

9.80* 

•H  W 
-P 

0 

>25 

0 

0 

23.20 

10 . 40 * 

* 

cc  •» 
•H  >> 

C 

>25 

0 

0 

24.70 

13.40. 

¥ 

2  u 

u  ca 

m 

•H  • 

>25 

0 

0 

26.20 

16.40. 

rH  U) 

CO  to 

e 

(U  0 

>25 

0 

0 

26.60 

17 . 20* 

U  0) 

o 

U  -H 

0)  -P 

>25 

0 

0 

27.60 

19 . 20 . 

* 

0  <D 
25  c 

,C  0 

•p  s 

>25 

0 

0 

28.40 

20.80 

>25 

0 

0 

29.10 

22.20 

250.0 

0.125 

31.20 

0 

0 

30.30 

24.60 

205.0 

0.145 

29.80 

0 

0 

32.00 

28.00 

158.0 

0.170 

26.80 

0 

0 

33.20 

30.40 

139.0 

0.195 

27.00 

0  4 

0 

34 . 70 

33.40 

111.0 

0.220 

24.40 

4.5x10  * 

0.0246 

35.30 

36.60 

92.80 

0.250 

23.20 

l.lxl0~3 

0.06 

37.80 

39.6 

79.50 

0.275 

21.80 

2.21x10  3 

0.12 

39.30 

42.60 

63.00 

0.325 

20.60 

3.5x10“;: 

0.190 

40.70 

45.40 

60.00 

0.340 

20.4 

3.7xlQ~3 

4xl0~3 

0.200 

>41.2 

46.40 

58.60 

0.345 

20.00 

0.218 

♦Those 

values  were 

used  to 

compare 

with  calculations  using  Bagnold  formula 

as  shown  in  reference  (7) . 


Table  9  (Con't) 


North 

Northwest 

Wes  t 

Southwest 

Q  total 

£  =900’ 

£  =  1700’ 

£  =  1400 

£  =  1400' 

inland 

r 

r 

r 

r 

t  hrs  Q 

(8)  (9) 

t  hrs  Q 

(8)  (9) 

t  hrs  Q 

(8)  (9) 

t  hrs  Q 

(8)  (9) 

in  lb  year 

0) 

u 

0) 


<D 

JS 

U 

p 

TO 

a> 

W 

o 

c 

c 

o 

■H 

• 

•H 

tn 

c 

4J 

0 

0} 

•H 

rH 

>> 

p 

3 

0 

U 

C3 

s 

r—i 

(Q 

a 

W 

0 

u 

a> 

C 

o 

0 

a> 

(0 

2 

c 

•H 

8 

4.4xl0"4 

10 

9.55X104 

1 

2 

2.7xl04 

14 

35. 5xl0_ 

1.52x10*? 

5 

0 

8 

2.03xlQ5 

3 

0 

8 

3.4xl05 

5 

4.00x10 

0 

3 

1.35x10^ 

2 

1.69X105 

9.2xl04 

0 

5 

2.44x10“ 

1 

0 

0.86xl04 

1 

0.31xl04 

15.12x10' 

0 

3 

2.24xl04 

40.44x10 

0 

0 

o  4 

35.50x10 

0 

3 

7.1x10 

81.10x10 

0 

0 

0 

30.40x10 

0 

0 

0 

33.6x10 

Total  = 


Z  2.3616xl06 
#/yr 


